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' Results are reported from the complete salt phase of the Sudbury Neutrino Observatory experiment in which 

NaCl was dissolved in the D2O target. The addition of salt enhanced the signal from neutron capture, as com- 
pared to the pure D 2 detector. By making a statistical separation of charged-current events from other types 
based on event-isotropy criteria, the effective electron recoil energy spectrum has been extracted. In units of 10 6 
cm' 2 s" 1 , the total flux of active-flavor neutrinos from 8 B decay in the Sun is found to be 4.94^ 2 |(stat)^ ^(syst) 
and the integral flux of electron neutrinos for an undistorted 8 B spectrum is 1 .68^JJ "^(stat)^" Qg(syst); the signal 
from (v v ,e) elastic scattering is equivalent to an electron-neutrino flux of 2.35^22(stat)^ 'jj(syst). These results 
are consistent with those expected for neutrino oscillations with the so-called Large Mixing Angle parameters, 
and also with an undistorted spectrum. A search for matter-enhancement effects in the Earth through a possible 
day-night asymmetry in the charged-current integral rate is consistent with no asymmetry. Including results 
from other experiments, the best-fit values for two-neutrino mixing parameters are Am 2 = (8.0^ *) x 10~ 5 eV 2 
and 6 = 33. 9^ 24 degrees. 

PACS numbers: 26.65. +t, 14.60.Pq, 13.15. +g, 95.85.Ry 



I. INTRODUCTION 

Results from the completed second phase of the Sudbury 
Neutrino Observatory (SNO) [1] are presented in this paper. 
The second phase began in June of 2001 with the addition of 



-2000 kg of NaCl [2] to SNO's -1000 tonnes of D z O, and 
ended in October 2003 when the NaCl was removed. The ad- 
dition of the salt enhanced SNO's ability to detect solar 8 B 
neutrinos in three ways. First, the neutron capture efficiency 
increased by nearly three-fold, allowing a statistically pre- 
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cise measurement of the neutral-current (NC) disintegration 
of deuterons by solar neutrinos. Second, the total energy of 
the y-rays from the neutron capture on 35 C1 is 2.32 Me V above 
the energy of the single y from the capture on deuterons. This 
higher y-ray energy approximately corresponds to a 1 MeV 
upward shift in the observed energy peak for neutrons and al- 
lows a precise measurement that is well above the low energy 
radioactive backgrounds to be made. Lastly, and perhaps most 
importantly, the isotropy of the photomultiplier tube (PMT) 
hit distribution on the geodesic array from multiple y-rays 
emitted after neutron capture on 35 C1 is significantly different 
from that produced by Cherenkov light emitted by a single 
relativistic electron. Hence, neutrons from the NC reaction 
and electrons from the charged-current (CC) interaction can 
be separated statistically without any assumptions about the 
underlying neutrino energy spectrum. 

Results from the first SNO phase [3-5], using pure D2O in 
the target volume, confirmed earlier measurements [6-9] of an 
observed deficit of solar electron-type neutrino flux compared 
to solar model expectations [10, 11] but additionally demon- 
strated, through measurement of the total active 8 B solar flux, 
that neutrinos undergo flavor transformation in transit to ter- 
restrial detectors. 

Initial measurements of the total active 8 B solar flux, based 
on the first 254.2 live days of the salt data set, have been 
published [12] and confirm and improve on results from the 
first D2O phase measurements [3-5]. In particular, the sta- 
tistical discrimination of CC and NC events with salt enabled 
an independent measure of the total active 8 B flux. The mea- 
sured flux is in very good agreement with solar model calcu- 
lations [10, 11, 13, 14] . 

The favored interpretation of these results is that neutri- 
nos undergo oscillations between flavor states given by lin- 
ear combinations of the non-degenerate mass eigenstates as 
described in the Maki-Nakagawa-Sakata-Pontecorvo (MNSP) 
framework [15]. For the case of 8 B solar neutrinos, the 
measurements support the Mikheyev-Smirnov-Wolfenstein 
(MSW) [16] hypothesis of matter enhanced oscillation, where 
electron neutrinos experience an additional interaction, com- 
pared with muon or tau neutrinos, in the presence of matter 
that can enhance neutrino oscillations. The SNO data, when 
combined with other solar neutrino measurements and reac- 
tor antineutrino results from the KamLAND experiment [17], 
show that neutrino oscillations are the dominant cause of fla- 
vor transformation and significantly restrict the allowed range 
of the relevant neutrino mixing parameters. 

In terms of neutrino mass and mixing parameters Am 2 and 
tan 2 9, solar neutrino data favor the so-called Large Mixing 
Angle (LMA) region. Maximal mixing is ruled out with a 
high-degree of confidence. 

The present paper extends the analysis to a total of 391 live 
days of SNO data from the salt phase, provides new results 
for the integral fluxes, and provides the CC energy spectrum 
and day-night spectral asymmetries. The values for tan 2 6 
and Am 2 are updated using a two-neutrino oscillation anal- 
ysis. The paper also provides a more detailed description of 
the full data analysis process for the SNO salt phase. 

The layout of the paper is as follows: Sec. II of the pa- 



per describes details of the detector hardware and software 
simulation. The data set, live time determination and event 
selection are discussed in Sec. III. 

The following sections discuss detector response and back- 
grounds, along with their associated systematic uncertainties 
on the neutrino measurements. Detailed systematic uncer- 
tainty evaluations are presented in the discussions of optical 
and energy calibration in Sec. IV, event vertex and direction 
reconstruction and isotropy in Sec. V, and neutron response 
in Sec. VI. Complete analyses of the many potential back- 
ground sources are given in Sec. VII followed by the evalu- 
ation of specific systematic uncertainties associated with the 
day-night asymmetry measurement in Sec. VIII. 

The procedure for analyzing the solar neutrino signal is 
discussed in Sec. IX. Solar neutrino results are presented 
with particular emphasis on the CC energy spectrum with the 
evaluation of differential energy systematic uncertainties in 
Sec. X, integral flux in Sec. XI, and day-night asymmetry in 
Sec. XII. Interpretation of the results in the context of the 
MSW framework for SNO data only and for the combined so- 
lar and reactor analysis can be found in Sec. XIII. A summary 
provided in Sec. XIV concludes the paper. 



II. SNO DETECTOR AND SIMULATION 
A. Detector 

The Sudbury Neutrino Observatory is a real time heavy wa- 
ter Cherenkov detector located in the Inco Ltd., Creighton 
mine near Sudbury, Ontario, Canada. The center of the de- 
tector is at a depth of 2092 m, or 6010 meters of water equiva- 
lent. At this depth, approximately 65 muons enter the detector 
per day. The neutrino target is 1000 tonnes of 99.92% isotopi- 
cally pure D2O contained inside a 12-m diameter acrylic ves- 
sel (AV). An array of 9456 20-cm Hamamatsu R1408 PMTs, 
which is mounted on an 18-m diameter stainless steel geodesic 
structure, is used to detect Cherenkov radiation in the target. 
A non-imaging light concentrator is mounted on each PMT to 
increase the effective photocathode coverage by the complete 
array to approximately 54% of An. 

To minimize the effects of radioactive backgrounds on the 
detection of solar neutrinos, materials with low intrinsic ra- 
dioactivity were selected for the construction of the detector. 
The acrylic vessel and the geodesic sphere are immersed in 
ultra-pure H2O to provide shielding against radioactive back- 
grounds from the geodesic structure and the cavity rock. An 
additional 91 PMTs are mounted looking outwards on the 
geodesic sphere and 23 PMTs are suspended facing inwards 
in the outer H2O volume to act as cosmic veto counters. Four 
PMTs that are installed in the neck region of the acrylic vessel 
provide veto signals to reject certain classes of instrumental 
background events (Sec. VII). Further details of the detector 
can be found in [1]. Note that for analysis purposes Cartesian 
coordinates are defined such that the center of the vessel is at 
(x,y,z) = (0,0,0), and the neck region is located symmetrically 
about the positive z axis. 

The SNO experiment detects solar neutrinos through the 
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charged-current (CC) and neutral-current (NC) interactions on 
the deuteron, and by elastic scattering (ES) on electrons: 

CC: v e + d -> p + p + e- - 1.442 MeV 
NC: v x + d -» p + n + v x - 2.224 MeV 
ES: v x + e~ — > v x + e~ 

where v x refers to any active flavor of neutrinos. The NC chan- 
nel has equal sensitivity to all active neutrinos, while the ES 
channel is sensitive primarily to electron-neutrinos. Hence, 
the NC measurement can determine the total active solar neu- 
trino flux even if electron-neutrinos transform to another ac- 
tive flavor [18]. 

In the first phase of the experiment with pure D2O, NC in- 
teractions were observed by detecting the 6.25-Me V y-ray fol- 
lowing capture of the neutron by the deuteron. For the second 
phase of data taking, (0.196±0.002)% by weight of purified 
NaCl was added to the D2O in May 2001 to increase the cap- 
ture and the detection efficiencies of the NC neutron. The ther- 
mal neutron capture cross section of 35 C1 is 44 b, which is sig- 
nificantly higher than that of the deuteron at 0.5 mb. When a 
neutron captures on 35 C1, the total energy released is 8.6 MeV. 
The combination of the increased cross section and the higher 
energy released results in a larger neutron detection efficiency 
at the same analysis threshold. 

Neutron capture on 35 C1 typically produces multiple y-rays 
(~2.5 per capture), while the CC and ES reactions produce 
single electrons. Each y-ray predominantly interacts through 
Compton scattering, producing an energetic electron. The 
Cherenkov light from neutron capture events, compared to 
that from CC and ES events, is more isotropic as the light is 
typically from several electrons rather than one. This greater 
isotropy, together with the strong directionality of ES events, 
allows good statistical separation of the event types. 

A precise measurement of the total active solar neutrino 
flux can be made through the NC channel without assump- 
tions about the underlying neutrino energy spectrum. This 
is relevant because the neutrino energy spectrum can be dis- 
torted from the generated 8 B spectrum via oscillation effects. 



B. SNO Monte Carlo Simulation 

The SNO Monte Carlo and analysis (SNOMAN) code is 
used for off-line analysis of the SNO data and provides an 
accurate model of the detector for simulating neutrino and 
background events. The Monte Carlo (MC) processor in 
SNOMAN provides processors for the generation of differ- 
ent classes of events, propagation of the primary particles and 
any secondary particles (such as Compton electrons) that are 
created, detection of the signal by the PMTs and simulation of 
the electronics response. With the exception of a few physics 
simulations (such as optical photon propagation), widely used 
packages such as EGS4 [19], MCNP [20] and FLUKA [21] 
are used in SNOMAN to provide accurate propagation of elec- 
tromagnetic showers, neutrons, and hadrons. 

Detailed models of all the detector components and cal- 
ibration sources are implemented in SNOMAN. Generators 



TABLE I: Tabulation of detector operational status. 



Data 


Percentage of 




run period 


Detector maintenance, 




v runs rejected by data quality checks, 




or elevated Rn or 24 Na levels 


22.6 


Calibration activities 


20.2 


Detector off 


6.2 


Selected v runs for analysis 


51.0 



for neutrino and calibration source signals, radioactive back- 
grounds and cosmic rays are also provided. Input parameters 
such as optical attenuation coefficients are determined from 
detector calibration. Calibration and detector parameters are 
input to SNOMAN and probability density functions (PDFs), 
used in the neutrino analysis, are generated. These features 
allow a direct assessment of the systematic uncertainties in 
physics measurements by comparing the detector responses 
for various calibration sources with the predictions of SNO- 
MAN. 

For the analysis of SNO data, SNOMAN provides various 
processors to unpack the data, to provide charge and time cal- 
ibration of the PMT hits for each event, to reconstruct event 
position and direction, and to estimate the event energy. 



III. DATA SET AND EVENT SELECTION 
A. Data Set and Live Time 

The measurements reported here are based on analysis of 
391.432 ± 0.082 live days of data recorded between July 26, 
2001 and August 28, 2003. As described below, 176.51 1 days 
of the live time were recorded during the day and 214.921 
days during the night. 

The selection of solar neutrino data runs for analysis is 
based upon the evaluation of detector operation logs and out- 
puts from an automated off-line data evaluation processor in 
SNOMAN. The automated processor checks the validity of 
event times, the trigger thresholds, the status of the electronics 
channels, and other detector parameters to ensure the quality 
of the runs selected for analysis. In addition, neutrino runs that 
are known to have elevated levels of radioactive background 
were removed from the data set. For instance, some neutrino 
runs have an elevated level of Rn ingress while others contain 
residual 24 Na radioactivity resulting from neutron activation 
during detector calibration. 

The bulk of the time that the detector was not live for neu- 
trino data acquisition was used for detector calibration and 
maintenance activities. Table I provides a tabulation of the 
detector operational status. 

The raw live time of the data set is calculated using a GPS 
synched 10 MHz clock on a run-by-run basis from the differ- 
ence in times between the first and last triggered events in each 
run. These calculated time intervals are verified by comparing 
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the results against those measured independently with the 50 
MHz detector system clock. Separate day and night live times 
are determined by splitting each run based on solar zenith an- 
gle Z where day = (cos 6 Z > 0) and night = (cos 8 Z < 0). The 
ratio of day to night live times is 0.82128 with an uncertainty 
of less than 5 x 10~ 7 . 

Several data selection cuts remove small periods of time 
from the data set. These include time intervals following 
high-energy cosmic-ray events and intervals containing time- 
correlated instrumental events To calculate the final live time 
for the neutrino data set, the total time removed by the full set 
of data selection cuts is subtracted. Data selection cuts remove 
a combined total of 1.8% of the raw live time for the neutrino 
data set. The dominant effect is caused by the cosmic ray veto 
cut. 

The live time calculation, including the corrections due to 
data selection cuts, is checked with an analysis of data from 
the detector diagnostic trigger. This pulsed global trigger 
(PGT) is a detector-wide trigger issued at a frequency of 5Hz 
based on timing from the 50 MH clock. Systematic uncer- 
tainties in live time are evaluated by comparing the PGT mea- 
surement to the 10MHz clock measurement, and by analyzing 
electronics and data acquisition effects that could prevent the 
detector from being live to neutrino data for short times. The 
total live time uncertainty is calculated to be ±0.021%. 



B. Event Selection 

The first step in selecting the solar neutrino candidate 
events involves the rejection of instrumental backgrounds and 
residual backgrounds from cosmic rays. A typical data rate 
was ~20 Hz, which was dominated by low-energy radioactive 
backgrounds in the detector and the PGT. The primary con- 
tributors to the instrumental backgrounds include events gen- 
erated by static discharge inside the PMTs, known as flashers, 
events produced when light is emitted from the neck region 
of the acrylic vessel, known as neck events, and electronic 
pickup events. The typical combined rate for these events 
is approximately one per minute compared to the ~10 per 
day rate of solar neutrino events. These instrumental back- 
ground events are identified and removed based on analysis 
of the charge and timing distributions of the triggered PMTs, 
event geometry in detector and in electronics space, signals 
from the cosmic veto counter outward-looking PMTs, PMTs 
installed to identify neck events, and time between events. 
Additional cuts using reconstruction information, referred to 
here as "high-level" cuts (Sec. VII), are used to remove events 
that do not possess the timing and isotropy characteristics of 
Cherenkov light from either single or multiple /?s. 

After the instrumental background and high-level cuts, two 
cuts are applied to remove cosmic -ray events. The first cut re- 
moves events that occur in a 20 s time interval after each event 
identified as a muon. The principal feature of muon identifi- 
cation is the requirement that there be at least 5 PMT hits in 
the outward-looking PMTs and 150 PMT hits (equivalent to 
~20 MeV) in the inner detector. A second, simpler cut re- 
moves any event following within 250 ms of any event with 
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FIG. 1 : Reduction of the data set as successive cuts are applied. The 
PMT instrumental, external light and pickup cuts remove instrumen- 
tal backgrounds originating from the detector hardware. The high- 
level cuts (Sec. VII) further reduce the instrumental backgrounds by 
rejecting events that do not possess the characteristics of Cherenkov 
light emission from single or multiple /3s. The fiducial volume cut, 
which selects events reconstructed with p < 0.77, removes most of 
the radioactive background events that originate outside the D2O tar- 
get. 



at least 150 PMT hits. This cut removes neutrons from muons 
missed by muon identification and neutrons induced by most 
atmospheric neutrino interactions inside the detector. 

The primary background sources are low in energy and 
external to the D2O volume. Hence, to minimize the sys- 
tematic uncertainties associated with backgrounds, all events 
selected for the results described here are required to have 
a reconstructed vertex position within 550 cm of the cen- 
ter of the detector. In terms of the volume-weighted vari- 
able p = (R/Rav) 3 , where /?av=600.5 cm is the radius of the 
acrylic vessel that contains the D2O, the fiducial volume cut 
is p < 0.77. Candidate events are required to have an effective 
electron kinetic energy T e s (see Sec. IV) greater than 5.5 MeV. 
Details of all background sources are provided in Sec. VII. 

The results of sequentially applying the cuts to the data set 
are shown in Fig. 1. The instrumental cuts are applied in steps 
to show the removal of various classes of backgrounds. The 
application of the high-level cuts then reduces the data set 
further, leaving the set of neutrino candidates. Some of the 
cuts remove individual events based upon their characteris- 
tics, while others remove periods of time from the data set. 
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TABLE II: Signal loss of the instrumental and high-level cuts for 
each signal class. These measurements are averaged over the energy 
spectra above r eff = 5.5 MeV. For CC and ES signals an undistorted 
8 B neutrino spectrum has been assumed. 



Signal Class 


Signal Loss, % 


Charged current 
Neutrons 
Elastic scattering 


0.57!™ 
0.68*»f 
0.86*$ 



The efficiency of the event-based cuts is measured inde- 
pendently for electron and neutron events. For neutrons, a 
software tagging approach (Sec. VI) is used to select neu- 
tron events from a 252 Cf fission source. These tagged neutron 
events are then used to measure the efficiency of all cuts. For 
electrons, the efficiencies of the instrumental cuts and high- 
level cuts are measured independently. The efficiency of the 
instrumental cuts is measured using tagged /? events from a 
8 Li source [22], while 6.13-MeV y-rays from a tagged 16 N 
source [23] are used to measure systematic uncertainties. The 
efficiency of the high-level cuts is established by MC simula- 
tions. Tagged 16 N events are again used to measure systematic 
uncertainties and corrections. The estimated signal loss for 
each class of event, integrated over the expected distributions, 
is shown in Table II. 

Signal loss as a function of energy for CC electrons is 
shown in Fig. 2. This figure shows separately the system- 
atic uncertainties that are correlated between energy bins and 
those that are independent. 

Neutrons produced by cosmic-ray muons and atmospheric 
neutrinos are used to measure the day-night asymmetry in the 
cut efficiency. Because muon-induced neutrons are generated 
throughout the live time of the experiment, they correctly sam- 
ple any variations in cut efficiency as a function of time. Can- 
didate neutron events are identified by selecting events with 
a reconstructed vertex radius less than 550 cm, energy in the 
range from 6 to 10 MeV, and inside a time window between 
4 /is to 40 ms after any muon. Note that the mean neutron cap- 
ture time is 5.3 ms at the center of the detector and essentially 
all neutrons are captured by 40 ms. 

The instrumental background cuts are applied to produce a 
clean sample of muon-induced neutrons, and the efficiency of 
the high-level cuts is estimated for the sample. Alternatively, 
the high-level cuts are applied and the cut efficiency of the 
instrumental background cuts on the remaining events is esti- 
mated. The loss of neutron events from all cuts, as measured 
with muon-induced neutrons, is 0.94 ±0.1 7% and is consistent 
with the results from neutron calibration sources. The mea- 
sured day-night asymmetry, defined as the difference between 
the night and the day signal loss divided by their average, is 
0.18 ±0.33. 

After all data selection cuts have been applied, including 
the energy threshold and fiducial volume cuts, 4722 candidate 
neutrino events remain. 




r eff (MeV) 

FIG. 2: Signal loss as a function of energy for CC electrons assuming 
an undistorted 8 B spectrum. Signal loss uncertainties are divided into 
two classes: uncorrelated (error bars) and correlated (error bands). 
Correlated uncertainties arise from systematic uncertainties in the 
measurement of signal loss and uncorrelated uncertainties arise from 
statistical uncertainties in the calibration data used in this measure- 
ment. 



IV. DETECTOR CALIBRATION 

Interpretation of SNO's signals requires measurement and 
calibration of the detector components and response. Many of 
the details of the calibration of the detector components are 
described in [1] and will not be discussed here. 

A variety of calibration sources are deployed in the heavy 
and light water regions to characterize the detector response. 
Source deployment is achieved with a manipulator system [1] 
that is able to maneuver sources to various positions in two 
orthogonal planes within the D2O and along 6 vertical lines in 
the H2O. The positional accuracy of the manipulator system 
is ~2 cm along the central axis and ~5 cm off-axis in the D 2 0. 
In the H2O the accuracy is ~2 cm. 

Table III lists the primary calibration sources used. These 
include pulsed nitrogen laser light for optical calibration and 
PMT timing, 16 N y-rays to produce a reliable energy cal- 
ibration, 8 Li for energy and reconstruction calibration, a 
3 H(p, y) 4 He source ("pT source") to test linearity of the en- 
ergy scale, 252 Cf and Am-Be sources of neutrons, and U, Th, 
Rn, neutron activated 24 Na, and 88 Y to test detector response 
to backgrounds. Note that the "Rn spike" was a controlled 
release of a measured quantity of Rn gas into the D2O. 

The following section describes the optical and energy cal- 
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ibrations and our evaluation of the systematic uncertainties 
associated with these assessments. Later sections provide de- 
tails on event reconstruction, and determination of the detector 
response to neutrons and backgrounds. 



TABLE III: Primary calibration sources. 



Calibration source 


Details 


Calibration 


Ref. 


Pulsed nitrogen laser 


337, 369, 385, 


Optical & 


[1] 




420, 505, 619 nm 


timing calibration 




L6 N 


6.13-MeV y-rays 


Energy & reconstruction 


[23] 


8 Li 


/3 spectrum 


Energy & reconstruction 


[22] 


252 Cf 


neutrons 


Neutron response 


[1] 


Am-Be 


neutrons 


Neutron response 




3 H(p,y) 4 He ("pT") 


19.8-MeV y-rays 


Energy linearity 


[24] 


U,Th 


P-y 


Backgrounds 


[1] 


88y 


P-y 


Backgrounds 




Dissolved Rn spike 


P-y 


Backgrounds 




In-situ 24 Na activation 


P-y 


Backgrounds 





A. Optical Calibration 

The detector optical response parameters are required for 
MC simulations and for the energy reconstruction processor. 
The optical parameters are determined by analyzing the detec- 
tor response to photons generated with a pulsed nitrogen laser. 
Laser light is transmitted from the laser through optical fibers 
to a diffusing ball that can be positioned at various locations 
in the detector. 

A complete optical scan typically consists of measurements 
taken at approximately ~ 40 different positions and at six 
wavelengths (337, 369, 385, 420, 505, and 619 nm). This 
set of wavelengths spans SNO's detectable Cherenkov light 
spectrum. 

The laser plus diffuser ball system produces short (0.8 ns) 
pulses of light in the detector. For any given position of the 
diffuser ball, the PMT array measures a time distribution of 
light. The difference between the measured time of each PMT 
hit and the expected time-of-flight from the source position to 
PMT is called the "time-residual." The time-residual distribu- 
tion exhibits a large peak originating from "prompt", or un- 
scattered, light. The full width at half maximum of this peak, 
after all timing corrections have been applied, is ~ 1.8 ns. 
Smaller peaks occur between 10 and 100 ns later. These "late" 
peaks arise from light reflected from various surfaces within 
the detector and from late pulsing of the PMTs. The optical 
model analysis is restricted to a ±4 ns window centered on the 
prompt peak to reduce sensitivities to the complicated struc- 
ture of the late light distribution. Additional analyses with a 
±10 ns prompt peak window are conducted to verify the re- 
sults and evaluate sensitivities to the timing window cut. 

The optical model is used to predict the number of prompt 
PMT counts in each laser calibration run based on operational 
parameters of the laser source and the detector optical param- 
eters. The relevant operational parameters of the laser source 
are its light intensity and angular distribution. The optical pa- 
rameters of the detector are the D2O attenuation length, the 
combined acrylic and H2O attenuation length, and the rela- 



tive efficiency of the PMT-reflector assembly as a function of 
incidence angle. These optical parameters are extracted by 
fitting the data collected during a multi-position scan in the 
D2O using ax 2 minimization method. Note that the fits return 
the inverse of the attenuation lengths or "attenuation coeffi- 
cients". This technique is not sensitive to the separate AV and 
H2O attenuation coefficients, only the combined AV+H2O at- 
tenuation coefficients. The AV attenuation coefficients, given 
in Table IV, were obtained from ex- situ measurements [25]- 
[26] and the H2O attenuation coefficients are determined from 
the difference between the AV attenuation coefficients and the 
measured sum. 

TABLE IV: Acrylic vessel attenuation coefficients. 



Wavelength (nm) Attenuation coefficient (10 i cm ') 

337 56^4 

369 23.0 

385 12.2 

420 7.70 

505 7.09 

619 7.09 



The measured D2O and rLO+acrylic attenuation coeffi- 
cients include the effect of Rayleigh scattering which removes 
a fraction of the light from the prompt time window. Because 
the MC simulation must model both absorption and scattering, 
the scattering contribution is subtracted from the measured co- 
efficient and the resulting absorption coefficient is used as an 
input to the MC. The scattering coefficient is determined us- 
ing laserball data taken with a collimating mask over the ball 
and measuring the hit probabilities for PMTs outside of the 
angular acceptance of the collimated beam. 

Figure 3 shows sample D2O and H2O attenuation coef- 
ficients, measured at 369 nm and 420 nm respectively, for 
six scans taken during the salt phase. The H2O coefficients 
are constant within the accuracy of the measurements, while 
the D2O values exhibit a steady increase until late in the salt 
phase. For both MC simulation and energy reconstruction, 
the D2O attenuation coefficients are determined, based on the 
date at which the given run was taken, from a linear fit to the 
measured attenuation coefficients as a function of time. 

Chemical assays of the D2O indicate that the change in at- 
tenuation can be attributed to trace levels of contaminants. In 
particular, the possible presence of organic complexes in the 
salt phase, as well as measured increases in Mn contamina- 
tion, likely associated with the MnOx assays discussed be- 
low, are correlated with the increase in the attenuation coeffi- 
cients. As described below, the changing response was inde- 
pendently measured by the laser and the 16 N sources, and was 
corrected for in the data processing. Following desalination 
by reverse osmosis and water purification after completion of 
the salt phase, the attenuation levels returned to those mea- 
sured in June 2000 prior to salt deployment. 

In addition to the attenuation coefficients, the relative av- 
erage response of the PMT-reflector assemblies as a function 
of photon incidence angle is measured. A typical angular re- 
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FIG. 3: Attenuation coefficients of (a) D2O and (b) H2O as a func- 
tion of date for 369 and 420 nm pulsed laser scans. The data points 
are correlated by systematics that are common to each of the mea- 
surements. Note that the H2O values are determined by subtracting 
the acrylic vessel ex-situ measured attenuation coefficients from the 
measured H 2 0+acrylic values. 

sponse distribution, given in Fig. 4, presents angular response 
measurements for 385-nm light from three laser scans taken 
during the salt phase. The response is normalized to unity at 
normal incidence angle and, as can be seen, the response in- 
creases by ~ 12% at 35°. The angular response curves are 
reasonably constant over the salt running period and the aver- 
age of the response functions from all laser scans is used as 
input to the MC simulation and the energy reconstruction. 



B. Energy Calibration 

Once the optical calibration constants are determined, the 
energy calibration and systematic uncertainty evaluations are 
carried out. 

The simplest energy estimator in SNO is the number of 
PMTs that trigger in an event (/V m ts)- However, for a given 



FIG. 4: Measured angular response curves for PMT-reflector assem- 
blies from three optical scans taken at 385 nm during the salt phase 
of operation. The incident angle for photons originating from the 
fiducial volume (p < 0.77) is confined to less than 35 degrees. The 
average response from all the laser scans is used as the input to the 
MC simulation. Note that the y axis zero is suppressed. 

event energy, the corresponding mean number of triggered 
PMTs varies with event position and direction owing primar- 
ily to the effects of the D2O, AV, and H2O attenuation coeffi- 
cients and the varying PMT angular response. 

The late light, described above, is difficult to accurately 
model. In order to minimize systematic uncertainties asso- 
ciated with the late light, only the number of PMTs that fire 
within the "prompt" range of +10 ns centered on the time- 
residual peak is used to estimate the event energy. Figure 5, 
produced from 16 N events generated at the center of the detec- 
tor, illustrates the shape of a typical time-residual distribution 
for Cherenkov light. 

The number of prompt PMT hits (7Y pro mpt) is corrected for 
noise (/V no i se ), optical response (e respo nse) relative to response at 
the detector center (eo), and for the fraction of working PMTs 
at the time of the event (PMr wor ki ng ) compared to the total 
number of PMTs {PMT toVd \ — 9456) to produce the corrected 
variable TYcorrected: 



N, 



corrected — 



(^prompt ~ jVnoise) PMT lotz \ 

ie /e PMr working ' 



(1) 



c responsi 



This is the effective number of prompt PMT hits that would 
have fired in an ideal detector with the event vertex at the cen- 
ter of the D2O volume. During the salt phase of operation, 
the number of working PMTs was generally between 8800 
and 8600. The correction for PMT dark noise, measured typi- 
cally to be ~0.1 hits/event, is small when compared to the av- 
erage response for 5.0-MeV kinetic energy electrons of ~35 
corrected hits/event. 

The optical response function corrects for the relative ef- 
fects of path length through D2O, AV, and H2O and for inci- 
dence angle onto the PMT-reflector assembly compared to an 
event at the center of the detector. The response function is 
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FIG. 5: PMT time-residual spectrum for a 16 N run taken with the 
source positioned at the center of the detector. The time residual is 
calculated without any walk correction to the PMT hit time. Only hit 
PMTs with time residuals within the 20-ns "prompt window" cen- 
tered at the "prompt peak" are used to estimate the energy of an 
event. 



With the global efficiency set, electron MC events are then 
generated at a series of fixed energies to produce a "look-up" 
table that translates A^ corre cted into equivalent electron kinetic 
energy T eS . 



given by 

Response = £ £ £ ^f^C/, 0, ff.f, A) (2) 

W </>' A 

xM(r, 6', <j)')g{9', ^e-^e-^-e-^ 

where the sums are over 10 polar ( ff ) and 10 azimuthal (</>') 
angle bins relative to the reconstructed event vertex and direc- 
tion (6'-Q), and wavelengths A in a range (220-710 nm) that 
encompasses the wavelengths to which the detector is sensi- 
tive. €pmt(A) is the average wavelength response of the PMT- 
reflector assembly, P(r,6,6' ,<p' , A) represents the angular re- 
sponse function, g(6',(p') the Cherenkov light weighting dis- 
tribution, M(r, 6', <p') a correction for multiple photon hits in 
the PMTs, and are A dependent attenuation factors for 
the three media (1 = D 2 0, 2 = AV, 3 = H 2 0). 

The next step in the energy calibration is to translate 
AWected into an energy. A combination of data and MC sim- 
ulated high-rate (~ 200 Hz) 6. 13 -MeV y-ray events from the 
16 N source and MC simulated electrons are utilized for this 
purpose. 

First, MC simulated 16 N events are generated with the 
source at the center of the detector. The optical constants de- 
scribed above are used as inputs to the simulation, leaving an 
overall global PMT efficiency factor free which is tuned by 
matching the mean of the MC simulated A^ cone cted distribution 
to the corresponding distribution for 16 N data, as shown in 
Fig. 6. This sets the MC simulated A^ correcte d scale to that of 
the SNO detector. 



C. Systematic Uncertainties 

Among the most important systematic uncertainties on the 
measured fluxes and energy spectra are the detector energy 
scale and energy resolution uncertainties. Near the analysis 
energy threshold of T e s = 5.5 MeV, these uncertainties are ex- 
tracted from comparisons between 16 N data and MC simula- 
tions. 16 N data are recorded at regular intervals throughout the 
salt phase. This consists of approximately monthly deploy- 
ments at the center of the vessel and periodic scans throughout 
the x-z and the y-z planes in the D2O volume. Determinations 
of energy scale and energy resolution systematics are tested 
and extended with analyses of 8 Li source data taken during 
the salt phase and pT source measurements taken during the 
D2O phase. 

The factor PMT tota i/ 'PMT^o^g in equation 1 directly af- 
fects the energy scale and requires accurate identification of 
PMTs exhibiting normal working behavior. Normal PMT op- 
eration is assessed from analysis of PMT charge and timing 
distributions from high-statistics laser calibration runs. Anal- 
ysis of PMT occupancy distributions from 16 N data is used 
evaluate the average fraction of improperly functioning PMTs 
missed in analysis of the laser data. The estimated energy 
scale uncertainty associated with misidentification of improp- 
erly working PMTs is conservatively estimated to be 0.20% 

The energy scale response of the detector is sensitive to the 
electronic threshold and PMT gain. Cross-talk measurements 
and the shape of PMT charge distributions are both sensitive 
to such changes and are evaluated with 16 N calibration data. 
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Special 16 N runs were taken during the salt phase in which 
the PMT high voltages were varied, and separately with the 
data acquisition threshold settings varied, to simulate the ef- 
fects of gain and threshold changes on the charge and cross- 
talk distributions. Comparison to the standard set of 16 N data 
places limits on possible threshold and gain effects on the en- 
ergy scale of 0.20% and 0.40% respectively. 

Differences in the detector response as a function of event 
rate can directly affect comparisons between I6 N source runs 
(high-rate) and neutrino data (low-rate). 16 N runs, with the 
source event rate tuned to simulate the neutrino mode rate, are 
recorded at approximately monthly intervals. Rate effects are 
evaluated by comparing low-rate and high-rate 16 N data and 
are estimated to be less than 0.1%. Additional energy scale 
uncertainties associated with timing resolution are evaluated 
to be less than 0.1%. 

The data span approximately two years of detector oper- 
ation. To accurately extract the integral flux and day-night 
energy spectra, it is critical to evaluate and model the time- 
dependence of the detector state. As indicated above, the D2O 
attenuation coefficients increased slowly throughout most of 
the salt phase. 

For each 16 N data and MC simulated run, estimates of the 
mean and width of the iVcorreeted an d energy distributions are 
generated by fitting a Gaussian function to the central por- 
tion of the spectra. Figure 7 shows the distributions of mean 
^corrected and reconstructed energy as functions of date for 16 N 
runs with the source positioned at the center of the detector. 
The MC simulated runs have been generated with D2O atten- 
uation coefficients increasing as shown in Fig. 3. As is seen, 
the decreasing MC simulated response in Af corre cted matches 
the slope observed in the 16 N data. 

During the final period of the salt phase, no MnOx assay 
(Sec. VII) was taken and measurements indicate that the de- 
tector energy response also ceased to change. The mean en- 
ergy distributions shown in Fig. 7(b) have been corrected for 
the changing attenuation. The run-by-run differences between 
data and MC simulation mean and width estimates are taken 
as the temporal stability systematic uncertainties on energy 
scale and resolution respectively. 

The energy scale uncertainty arising from the temporal sta- 
bility evaluation is 0.15%. The energy resolution stability un- 
certainty is determined to be 1.8% and is dominated by an 
average offset between the data and the MC simulation. This 
offset is attributed to a combination of effects from electronic 
cross-talk and from tube-to-tube variation in PMT efficiencies 
not modeled in the MC simulation. 

Modeling the detector and the calibration sources involves 
a variety of simplifications and uncertainties. For the 16 N 
source, uncertainties in the 16 N decay branching ratios, mi- 
nor differences between pure D2O and the salt brine in the 
production and propagation of Cherenkov light, exact details 
in describing the source geometries, the finite step size in the 
EGS4 simulation, neglecting the minor velocity dependencies 
in the wavelength spectra of Cherenkov light, and the wave- 
length dependence of the index of refraction are estimated, in 
total, to contribute a 0.65% systematic uncertainty to the en- 
ergy scale. 
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FIG. 7: (a) Mean /V c01Tecled and (b) mean T cf j versus date for data and 
MC high-rate I6 N calibrations runs with the source at the center of the 
detector. Note that the full y-axis ranges are ~ 4% and ~ 2% of the 
average N comae a and T c ff values respectively. Error bars are statistical 
only, and the spread of the variation between data and Monte Carlo 
provides the measure of the energy scale uncertainty arising from 
temporal variations in detector response. 



Significant contributions to the energy scale uncertainty 
arise from evaluations of the radial response and detector 
asymmetries. Figure 8 shows data and MC simulation mean 
energy values and their ratio versus source radial position p. 
The volume-weighted mean difference between data and MC 
simulation of 0.45% is taken as the radial response energy 
scale uncertainty contribution. 

The detector asymmetry component means point-to-point 
non-uniformities in detector response due to asymmetric de- 
tector features such as support ropes and the neck of the 
AV. It is evaluated by determining the volume-weighted av- 
erage standard deviation in each of seven radial bins of the 
data/MC simulation ratio distribution and is calculated to 
be 0.59%. The corresponding energy resolution radial and de- 



10 




1.05 




FIG. 8: (a) Data and MC mean energy versus p distributions and (b) 
the run-by-run ratio of data to MC mean energy versus p are shown 
for 16 N calibration runs. The dashed vertical line corresponds to the 
fiducial volume cut at R = 550 cm. Points at the same value of p can 
have differing energy response in data or Monte Carlo due to local 
point-to-point non-uniformities in detector response. 



the 16 N source. The effects of cross-talk and noise on the en- 
ergy scale are estimated to be < 0.25%. The energy resolution 
uncertainty was determined from pT data to be 10% at 19.8 
MeV. Hence the energy resolution uncertainty is applied as 
the function 

r eff < 5 MeV : Ao- T = 3.4% (3) 
r eff > 5 MeV : Acr r = [3.4 + 0.478(r eff - 5)]%. (4) 

The energy response for electrons was characterized as a 
Gaussian function with resolution cr T = -0.131 +0.383 
0.0373 \T e , where T e is the true electron kinetic energy in 
MeV. Table V summarizes all contributions to energy scale 
and resolution systematic uncertainties. Energy scale contri- 
butions are added in quadrature giving a total uncertainty of 
1.15%. The radial and detector-asymmetry energy resolution 
components are added together in quadrature and then added 
linearly to the data-MC offset to produce a 3.4% total uncer- 
tainty. 



TABLE V: Summary of energy scale and resolution systematic un- 
certainties. 



Scale uncertainty 



Source 


Uncertainty 


Detector PMT status 


0.20% 


Electronics threshold 


0.20% 


Electronics gain 


0.40% 


Electronics rate effects 


0.10% 


Time calibration 


0.10% 


Time drift/stability: data-MC 


0.15% 


Radial distribution: data-MC 


0.45% 


Detector asymmetry 


0.59% 


16 N source modeling 


0.65% 


Cross-talk/pickup non-linearity 


0.25% 


Total 


1.15% 


Resolution uncertainty 


Source 


Uncertainty 


Central lb N runs: data-MC 


1.8% 


Detector asymmetry 


1.4% 


Radial dependence 


0.8% 


Total 


3.4% 



tector asymmetry systematic uncertainties are evaluated in an 
analogous fashion and found to be 1.4% and 0.78% respec- 
tively. 

These estimates of energy scale and resolution uncertainties 
are obtained near the detector threshold. The cross-talk and 
multiphoton effects of higher energy events are probed with 
high-intensity pulsed laser data. In addition, the pT source, 
which generates 19.8 MeV y-rays, enables a direct test of the 
higher energy scale systematic uncertainties. However, high- 
rate of neutrons emanating from the pT source precluded its 
deployment during the salt phase of the experiment. Compar- 
ison of the pure D2O and salt laser runs indicate no additional 
unmodeled effects and supports the application of the pT data 
to the salt data set. It is found that the energy scale uncertainty, 
evaluated from pT data, is not greater than that evaluated with 



V. EVENT RECONSTRUCTION 

For each event, PMT trigger times and positions are used 
to reconstruct the event vertex, direction, and isotropy. The 
following sections outline the algorithms and the determina- 
tion of systematic uncertainties associated with these recon- 
structed variables. 



A. Event Vertex and Direction 

Event positions and directions are reconstructed by analysis 
of the times and positions of triggered PMTs in each event. 
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The time-of-flight corrected PMT trigger time is used to define 
the PMT time-residual 



T] eh - t\ - fflt - 



l^fit - n\ 

«eff 



(5) 



where f, and r t are the trigger time and the position of the i' h 
PMT in the event, and ffl t an d ?&t are the fit time an d the re- 
constructed position of the event vertex. The effective photon 
velocity, u e s = 21.87 cm/ns, is the group velocity of the mean 
detected photon wavelength at 380 nm in D2O. 

For each event a likelihood is constructed from the PMT 
time-residual probabilities where MC simulated events are 
used to derive the reference PDF. The PDF is approximated 
as a constant for time-residual greater than 15 ns because the 
time peaks associated with reflected photons strongly depend 
on the event location. A time cut of +50 ns relative to the 
median PMT hit time is imposed to reduce the effects of re- 
flected photons and PMT dark noise. Seed vertices are chosen 
randomly from the detector volume and the negative log like- 
lihood function is minimized with respect to fe t and r^ t until a 
global minimum is found. 

Event direction is reconstructed independently after the 
best-fit vertex has been found. It is estimated based on the 
assumption that the events produce Cherenkov light emitted 
in the characteristic cone shaped pattern. A likelihood is con- 
structed based on the distribution of directions from the recon- 
structed event position to the triggered PMTs relative to the 
corresponding MC calculated distribution for Cherenkov light 
events. A "prompt" time cut (+10 ns) on the time-residual dis- 
tribution excludes most Rayleigh scattered and reflected pho- 
tons. The negative log likelihood is minimized to determine 
the best fit direction. For multiple electron events the recon- 
structed direction tends to be weighted towards the direction 
of the most energetic electron(s) in the event. 

Vertex reconstruction uncertainties are evaluated by com- 
paring average reconstructed event positions of 16 N calibra- 
tion data with 16 N MC simulations and by comparing the av- 
erage reconstructed source position to manipulator estimated 
position for 16 N calibration data. The manipulator source po- 
sition measurement is most accurate when operating in single 
axis mode along the z-axis of the detector. Figure 9 shows the 
difference between mean reconstructed x, y, and z positions 
and manipulator estimated source position as functions of p. 
These figures indicate that the x and y vertex reconstruction 
uncertainties are not more than 2 cm. A larger difference in z 
is apparent and the z vertex reconstruction uncertainty is taken 
to be 6 cm. 

In addition to coordinate shifts, fiducial volume uncertainty 
is also evaluated. Radial scaling bias could occur through 
reconstruction biases or timing calibration uncertainty. Re- 
constructed radial distributions of events from calibration and 
neutrino signal data near the AV are compared to MC simu- 
lations to evaluate this uncertainty. The radial uncertainty is 
estimated to be 1.0% of the radius (i.e., 5.5 cm at 550 cm). 

Angular response uncertainty is determined from analysis 
of 16 N data and 16 N MC simulated events. The y-rays pro- 
duced by the 16 N source travel an average of 30 cm in D2O 
before a Compton scatter occurs. The vector from the 16 N y- 
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FIG. 9: Difference between mean reconstructed (a) x, (b) y, and (c) z 
coordinates and the coordinates of the deployed 16 N source versus p. 
The 16 N source was deployed along the central (z) axis for the data 
shown in this figure. From these data, systematic uncertainties on x 
and y coordinate vertex reconstruction are evaluated to be 2 cm and 
the uncertainty on z coordinate vertex reconstruction is evaluated to 
be 6 cm. 
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FIG. 10: Data and MC distributions of the cosine of the angle be- 
tween generated event direction and reconstructed direction for a 16 N 
source located at the center of the detector are shown. The events are 
selected such that the reconstructed vertex is at least 120 cm away 
from the source position. 

ray source to the reconstructed event vertex provides a good 
estimate of the generated electron direction, since for events 
above T e g = 5.0 MeV the Cherenkov light is dominated by 
that from Compton electrons that are forward-peaked. Gam- 
mas that travel at least 120 cm from the source before scatter- 
ing are employed to ensure the estimate of the event direction 
is reliable. For each data and MC simulated 16 N run, the dis- 
tribution of the cosine of the angle 9 between the electron di- 
rection and the reconstructed direction is fit with the function 



R = N[eP- 



i(cosfl-l) 



cos 0-1) 



(6) 



where N is the overall normalization, /3s parametrizes the dis- 
tribution for electrons scattered only a small amount, /3m ac- 
counts for those scattered through large angles, and the ratio 
of these components is Shown in Fig. 10 are sample data 
and MC simulated distributions of cos 9 for an 16 N run at the 
center of the detector. The angular resolution systematic un- 
certainty, determined by comparing the average difference be- 
tween data and MC simulation fit parameters for all 16 N runs 
taken during the salt phase, is 16% of /3s ■ 



B. Event Isotropy 

In the context of SNO event analysis and signal extraction, 
isotropy refers to the uniformity of the distribution of trig- 
gered PMTs on the PMT array. 

For data collected during the SNO pure D2O phase, the CC 
and ES signal events produced a single primary electron while 
the NC events produced a single 6.25-MeV capture y-ray. For 
a NC event above T e ff = 5.0 MeV the Cherenkov light is 
dominated by that from a single forward-scattered Compton 
electron, so all signal event types in the pure D2O phase had 
similar isotropy distributions. With salt added to the SNO 



detector the characteristic response to neutrons is multipho- 
ton. Event isotropy is measured from the spatial distribution 
of triggered PMTs and is an effective signal separation tool in 
this circumstance. 

Several variables constructed to measure isotropy are found 
to have comparable separation power between the single elec- 
tron (CC and ES) and the neutron (NC) signals. The variable 
used, which could be simply parameterized and facilitate sys- 
tematic uncertainty evaluations, is f3u = j3\ + \f3\ where 



N(N - 1) 



(7) 



i=l /=«+! 



In this expression Pi is the Legendre polynomial of order I, 
9jj is the angle between triggered PMTs / and j relative to 
the reconstructed event vertex, and N is the total number of 
triggered PMTs in the event. 

Initial comparisons of the isotropy distributions for I6 N data 
and MC simulated I6 N showed that the mean value of /3\\ for 
16 N data was ~2.5% larger than for the 16 N MC simulation. 
This was caused primarily by approximations used in the de- 
scription of electron scattering in the MC simulation. Electron 
transport is handled within the MC simulation by EGS4 [19] 
in which elastic collisions with atomic nuclei are modeled us- 
ing Moliere's theory of multiple-scattering as formulated by 
Bethe [27]. This neglects the effect of the spin of the elec- 
trons which slightly reduces the amount of large-angle scatter- 
ing. The MC simulation therefore overestimates the amount 
of Cherenkov light emitted in the backward hemisphere, pro- 
ducing a light distribution which is slightly too isotropic, i.e., 
a jSi4 slightly too small. There are also other approximations 
in the treatment of multiple-scattering in EGS4, that affect the 
distribution, but the effect of spin is the most significant. 

Including the effect of spin multiplies the Rutherford cross 
section by a factor M{9,/3,Z), first calculated by Mott [28], 
given by [29]: 
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M(9,/3,Z) = 1 -f3 L sin z - + n/3— sin -(1 - sin -), (8) 
2 Tic 2 2 

where 9 is the single scattering angle, /3 is the electron speed 
divided by c, and Z is the charge of the scattering nucleus. 
In EGS4, an electron is propagated a step length x at which 
point its direction is changed by an angle (p. This angle is 
drawn from a probability density distribution for electrons 
multiple-scattered by a screened Coulomb potential without 
the Mott terms. A simple MC program was written to evalu- 
ate the probability density distributions with and without the 
Mott terms, and correction constants were generated giving 
the amount Atp to be subtracted from the EGS4 angle tp for a 
number of different step lengths and electron kinetic energies. 
For large angles the correction is approximately step length 
independent and close to that obtained when assuming single 
scattering dominates, while for small angles and longer step 
lengths it is smaller. 

This correction to the EGS4 multiple-scattering angles was 
parameterized as a function of energy using the average EGS4 
step length for a given energy. The fraction of electrons that 
scatter into the forward hemisphere after passing through a 
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FIG. 11: j8,4 isotropy distributions for 252 Cf data and MC, I6 N data 
and MC, and simulated CC events. 
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layer of water was compared to that obtained with the up- 
dated version of EGS4, EGSnrc [30] which includes the Mott 
terms as well as other improvements. For electrons with a ki- 
netic energy of 5.0 MeV passing through 1 cm of water, the 
percentage increase in the fraction that scatters into the for- 
ward hemisphere over that obtained with EGS4 is 0.9% when 
applying the average step length correction, and 1.2% using 
EGSnrc. The correction to EGS4 was tuned to give the same 
percentage as EGSnrc. 

Figure 1 1 shows distributions of fiu using data from 252 Cf 
and 16 N sources and from corresponding MC simulations. 
Also shown is a MC simulated distribution for CC events. The 
252 Cf and CC events have an imposed kinetic energy thresh- 
old of 5.5 MeV, while the I6 N events have kinetic energies 
between ~ 4 MeV and ~ 6 MeV. The 16 N source emits 6.13- 
MeV y-rays, which undergo Compton scattering and produce 
more than one lower energy electron. The CC electrons un- 
dergo relatively less multiple-scattering per unit path length 
than the Compton scattered electrons from the 16 N source, 
as the CC electrons have higher energy. The average flu for 
16 N events is therefore smaller (more isotropic) than for CC 
events. 

The /?i4 parameter is correlated with event energy and to a 
lesser extent with radius. Multi-dimensional PDFs are there- 
fore used in the analysis, as described in the signal extraction 
section. 

Systematic uncertainties for j3u are evaluated through data- 
MC comparisons of 252 Cf and 16 N calibration runs. For mono- 
energetic data, /3i4 can be well approximated by a Gaussian 
function. Calibration data for 252 Cf and 16 N and the corre- 
sponding MC simulated distributions are fit run-by-run. The 
difference in the extracted means and widths of the fits is used 
to characterize the uncertainties. Temporal stability is mea- 
sured with runs taken with the sources at the center of the de- 
tector. Radial uncertainty is determined from multi-axis scans. 
Figure 12 shows the distributions of f}\n width for 252 Cf data 
and MC as a function of p, and /3u mean for 16 N data and MC 
versus p. The average temporal and volume-weighted radial 
data-MC differences are evaluated separately and then added 
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FIG. 12: f} u isotropy distributions of (a) 252 Cf data and MC width 
versus p and (b) I6 N data and MC mean versus p. Systematic effects 
observable in the distributions are attributed to temporal and spatial 
non-uniformities in detector response and give the estimates on the 
/?14 systematic uncertainties. 

in quadrature. Table VI gives the estimated systematic mean 
and width f3\n uncertainties from both 252 Cf and 16 N sources. 
These uncertainties are propagated through the signal extrac- 
tion by shifting the means and smearing the fin PDFs. 



TABLE VI: Summary of /?i4 scale and resolution systematic uncer- 
tainties. 



Uncertainty 


Source 


Scale Resolution 


252 Cf 


0.48% 0.67% 


W N 


0.85% 0.94% 



The energy dependence of the systematic uncertainty on/^4 
for CC events was evaluated using 8 Li calibration data in ad- 
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FIG. 13: Neutron energy response during the D 2 (dashed) and salt 
(solid) running periods. The vertical line represents the analysis en- 
ergy threshold of T eff = 5.5 MeV in the salt period. For the D 2 
period the analysis energy threshold was T e s = 5.0 MeV. The distri- 
butions shown here are normalized to the neutron detection efficiency 
in the two phases for R< 550 cm. 

dition to 16 N. 



VI. NEUTRON RESPONSE 

With salt added to the D2O volume, neutron capture is dom- 
inated by capture on 35 C1 and losses due to capture on 'H 
and I7 are significantly reduced. For a 252 Cf source at the 
center of the detector, the probabilities of neutron capture are 
90%( 35 C1), 4% ( 2 H), 2.5%('H ), with the remaining 3.5% ab- 
sorbed by oxygen, sodium or other isotopes (only 0.3% cap- 
ture on 37 C1). In D2O with no salt additive, the capture prob- 
abilities are 49% ( 2 H) and 30%('H), with 14.5% absorbed by 
oxygen isotopes in the heavy water and the remaining 6.5% 
escaping the D2O volume. For neutrons generated uniformly 
in the heavy water, the probability of capture on 35 C1 in the 
salt phase is about three times larger than that of capture on 
deuterium in the absence of salt. 

Not only is the capture efficiency increased, but the en- 
ergy deposited in the detector is also increased. As shown 
in Fig. 13, the peak of the energy distribution moves to higher 
energy so that, for the same energy cut, the salt phase has 
improved neutron detection efficiency compared to the D2O 
phase. This allowed a higher energy threshold for the salt 
phase and hence less low-energy background contamination. 

Neutron response is calibrated primarily with neutrons pro- 
duced by a 252 Cf source with secondary checks made by anal- 
ysis of neutrons generated by an Am-Be source and by MC 
simulations. To determine the neutron detection efficiency us- 
ing the 252 Cf fission source, the absolute neutron production 
rate (source strength) has to be determined. Shown in Table 
VII are the results of four different techniques used to evalu- 
ate the source strength and presented as inferred strength of 
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FIG. 14: Neutron capture efficiency versus radial position of the 
252 Cf source for the pure D 2 and salt phase. The solid line is a 
fit of the salt phase data to Eq. (10), and the dotted line is a fit of the 
D 2 phase data to a neutron diffusion model. 

the source on June 12, 2001. The 252 Cf activity decays with a 
half-life of 2.645 years, and this decay, together with that of a 
small 250 Cf contamination, is taken into account in evaluating 
the source strength at the time of a given calibration run. 

The Frisch Grid and triggered Si(Li) methods both use an 
array of calibrated 3 He detectors to detect neutrons with the 
Si(Li) method being triggered on the fission daughter products 
of 252 Cf. These methods provide independent measures of the 
source strength prior to deployment in the SNO detector. 

In-situ measurements of the source strength and efficiency 
have also been made in the D2O and dissolved-salt phases. 
The D2O multiplicity method is an in-situ method used in the 
pure D2O phase of SNO operation to determine the detection 
efficiency and fission rate. In this method the distribution of 
the number of neutrons detected in 2-second time windows is 
plotted and then fit to the multiplicity function. For a model 
in which the neutron capture time is negligible compared to 
this time window, the probability of detecting d neutrons in a 
given time window for / neutrons generated is given by 



P{d) = e- XI 5 d $ 

l=d 



/! 



d\{l-d)\ 



€ d (\ - €)'- d 



-(/-AW 2 

Ze 2Nc- 2 
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(2nNcr 2 ) 1 ^ N\ 



-AT 



(9) 



The neutron detection efficiency e and the fission rate A are 
the free parameters in the fit. The factor 6d,o is 1 for d — 
and is otherwise. The time window is T. The multiplicity 
of the 252 Cf source, fx, is taken to be 3.7676±0.0047 neutrons 
per fission, and the width of the multiplicity distribution, cr, 
is 1.57 [31]. The corrections due to the finite neutron capture 
lifetime in the D2O volume were estimated by MC simula- 
tions. 

A time-series method was used in the salt phase to extract 
the neutron detection efficiency and the fission rate by using 
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the time separation between fission daughter y-rays and neu- 
trons in a 5 -parameter fit to an analytical distribution that gen- 
eralizes the model on which Eq. (9) is based to explicitly in- 
clude the neutron capture lifetime and the effect of detecting 
a small fraction of the fission y-rays that accompany neutron 
production. As can be seen in Table VII, the various tech- 
niques are in agreement and we calculate a weighted mean as 
our best estimate of the neutron source strength. 



TABLE VII: Results of the various methods for determining the 252 Cf 
source strength. The source strength is determined for June 12, 2001 . 
The^ 2 for the tabulated source strengths is 5.6. 



Method 


Source strength (neutrons per second) 


LANL Frisch Grid 


16.75 + 0.14 


LANL triggered Si(Li) 


17.08 ± 0.43 


D 2 multiplicity 


16.33 ±0.18 


Salt time-series 


16.46 ±0.12 


Weighted mean 


16.55 ± 0.08 



The neutron efficiency is determined by comparing the 
number of neutrons detected to the number produced, as a 
function of the position of the 252 Cf source and the energy 
threshold. In this comparison it is important to take into ac- 
count events from fission y-rays that are emitted in the spon- 
taneous decay of 252 Cf. First the small fraction due to y-rays 
above T e g = 6.5 MeV is determined; then the neutron energy 
spectrum in the salt phase (Fig. 1 3) is used to extrapolate down 
to zero threshold to give the capture efficiency (Fig. 14), and 
to r e ff = 5.5 MeV (Fig. 13) to give the detection efficiency. 

To extract the shape of the 252 Cf y spectrum above r ei f = 
6.5 MeV, a 252 Cf source run was recorded, during the short 
pure D 2 phase following the removal of salt. In D 2 0, neu- 
trons can be separated easily from the source y-rays by re- 
quiring that neutron candidate events reconstruct more than 
150 cm from the source. The fission y-ray energy distribu- 
tion is then obtained by subtracting the energy distribution of 
the selected neutron events from the energy distribution of the 
events that reconstruct within 150 cm of the source. The y- 
ray energy distribution obtained from this measurement and 
the neutron energy distribution obtained from the salt phase 
(see below) are then fit to the raw 252 Cf energy distribution in 
salt with the scaling on the neutron and y distributions left as 
free parameters to obtain the y background correction above 
r ei f = 6.5 MeV. The time series analysis also gives an inde- 
pendent estimate for the gamma fraction. The combined result 
from these two methods yields (1.34+Jjg)%. 

To determine the neutron energy spectrum in the salt phase, 
a clean neutron sample is required. As the capture distance for 
neutrons in the salt phase is similar to the attenuation length 
of fission y-rays, a radial cut cannot be used to select neutron 
events from the calibration data. Instead a "burst cut" was 
developed to select neutron events from calibration data us- 
ing the coincidence between fission y-rays and neutrons. In 
salt, the mean capture time for neutrons from 252 Cf data at 
the detector center is measured to be 5.29 + 0.05 ms and af- 
ter approximately 40 ms almost all neutrons have been cap- 
tured. The mean time between fission bursts for the 252 Cf 



source used in SNO is about 250 ms. In the burst cut a fis- 
sion y-ray candidate event is selected by choosing events with 
no preceding events within a 50-ms time interval, and events 
in a time interval of 40 ms after the selected first event are 
tagged as neutrons. The burst cut has a selection efficiency for 
neutrons of 40%, but less than 0.1% of the selected neutron 
candidate events are estimated to be fission y-rays above the 
threshold (T e g = 5.5 MeV) used in this analysis. 

The capture efficiency is obtained for each 252 Cf run with 
source radial position r (start radius of the neutrons). These 
point source efficiency results are fit to the empirical model 



e(r) = A {tanh [B (r - C)] - 1) 



(10) 



where e(r) is the neutron capture efficiency at source position 
r and requiring neutrons to be captured (reconstructed vertex) 
inside R< 550 cm. A, B, and C are the fit parameters of the 
model. The volume-weighted capture efficiency 1 is then ob- 
tained from the ratio of integrals 



e = 



(11) 



where e(r) is plotted in Fig. 14. 

Figure 15 shows the comparison between detection effi- 
ciency distributions obtained from calibration data and from 
the NC MC simulation as functions of p after the r e tf > 
5.5 MeV selection criterion has been applied. The neutron 
detection efficiency along with its uncertainty derived from 
252 Cf calibration data are shown as the shaded band in this fig- 
ure. The ratio of the neutral current MC simulation efficiency 
to the efficiency obtained from calibration data is within 2.2% 
of one; which is within the estimated systematic uncertainty. 

A list of corrections applied to the 252 Cf efficiency mea- 
surement is summarized in Table VIII. Apart from the gamma 
fraction correction discussed previously, these corrections are 
calculated through MC simulation studies The source sam- 
pling correction is needed due to detector asymmetry. The 
252 Cf calibration data were predominantly collected with the 
source positioned in the bottom half of the detector, and there- 
fore do not sample the whole detector. This correction is cal- 
culated by doing MC simulation studies at the same 

252 cf 

source positions as the data, and comparing the efficiency 
from the point source MC simulation studies with the effi- 
ciency derived from a MC simulation study of 252 Cf neutrons 
uniformly distributed in the detector. The uncertainty on this 
correction is included in Table IX. A correction is applied 
to the calibration efficiency measurement to account for the 
2 H(n,2n)'H reaction with fission energy neutrons. This cor- 
rection is determined from a MC simulation study. The y 
fraction correction is the previously discussed correction to 
account for 252 Cf fission y-rays contaminating the 252 Cf neu- 
tron data. A correction is also applied to account for the 
16 0(n,a:) 13 C reaction with fission energy neutrons. Finally, 
a source geometry correction is applied to account for neutron 
captures on the steel and acrylic of the source holder. 

The systematic uncertainties on the 252 Cf neutron detection 
efficiency measurement are listed in Table IX. Some of the 
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TABLE IX: Systematic uncertainties on the neutron efficiency mea- 
surement from the calibration data. 



Source 



Uncertainty, 



Source strength ±0.5 

Source position +1.7, —1.0 

Gamma fraction +0.56, -1.05 

AV position ±0.3 

(n,2n) ±0.10 

(n,a) ±0.13 

Empirical fit - polynomial fit +0.4 

Source sampling ±1.0 

Source geometry ±0.53 



Total uncertainty 



+2.3, -2.0 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



FIG. 15: Comparison of neutron detection efficiency (T e s > 
5.5 MeV) for MC simulated NC events (data points) and that derived 
from 252 Cf calibration data (shaded band) as a function of volume- 
weighted radius p. The band represents the statistical and systematic 
uncertainties summarized in Table IX. An additional 1.0% radial re- 
construction uncertainty that is assigned to the solar neutrino flux is 
also included in the band. The volume-weighted NC MC efficiency is 
within the systematic uncertainty assigned to the 252 Cf measurement. 



uncertainties arise from the corrections listed in Table VIII. 
The source strength uncertainty is derived from the results 
summarized in Table VII. Source position uncertainty is ob- 
tained by shifting the estimated source positions by ±2 cm or 
±10 cm in radius, depending on the detector region, and then 
re-calculating the volume-weighted efficiency. The AV posi- 
tion uncertainty is taken from a +6 cm z shift in the acrylic 
vessel position. An estimate of the uncertainty in the interpo- 
lation shown in Fig. 14 was taken as the difference between a 
high-order polynomial fit and the empirical fit to the data. The 
volume-weighted neutron detection efficiency for the analy- 
sis threshold of T e s = 5.5 MeV and a fiducial volume of 
550 cm after applying the corrections listed above is (40.7 
± 0.5 ^o'g)% where the first uncertainty is statistical and the 
second the combined systematic uncertainty. This is the effi- 
ciency for detecting neutrons uniformly generated within the 
whole AV and with reconstructed vertex less than p - 0.77. 



TABLE VIII: Corrections applied to the neutron efficiency measure- 
ment from the calibration data. 



Source 


Correction, % 


Source sampling 


-(2.4 ± 1.0) 


(n,2n) 


-(0.58 ±0.10) 


(n,a) 


+(0.66 ±0.13) 


Gamma fraction 
Source geometry 


-(1.34!?;*) 
+(2.03±0.53) 


Total 


-1 73 +1 1 



VII. BACKGROUNDS 

Several sources of backgrounds are present in the data. 
These include instrumental backgrounds, backgrounds from 
the natural 232 Th and 238 U radioactivity chains, and back- 
grounds associated with products of cosmic muon spallation 
and atmospheric neutrino interactions in the detector. While 
some of these background types can be eliminated by analy- 
sis cuts, most cannot be distinguished from the solar neutrino 
signals. Table X provides a summary of the estimated contri- 
butions from these backgrounds. In the following subsections, 
the identification and the determination of the contributions of 
these backgrounds is discussed. 



A. Instrumental Backgrounds 

A significant portion of the events comprising the raw data 
are the instrumental backgrounds discussed in Sec. IIIB. The 
instrumental background cuts and the high-level cuts are very 
efficient at removing these events. The residual contamina- 
tion of instrumental backgrounds in the data set is measured 
by using a bifurcated analysis [32]. In this analysis, each set 
of cuts is used to calibrate the acceptance for background of 
the other set, allowing the leakage through the combination 
of both sets to be calculated. The two sets of cuts must be 
orthogonal (uncorrected) for the bifurcated analysis to work 
and are chosen appropriately. Orthogonality is demonstrated 
using a technique known as "relaxing the box" in which the 
bifurcated analysis correctly estimated the increase in residual 
background as cuts are relaxed when applied to the data set. 

The bifurcated analysis provides an upper limit on the resid- 
ual instrumental contamination of 3.0 events in the neutrino 
data which is treated as a 68% CL limit in subsequent anal- 
yses. This analysis is represented in Fig. 16 where the dis- 
tribution of the high-level cut parameters is presented for the 
neutrino data set before application of these cuts, and for in- 
strumental backgrounds rejected by the instrumental cuts. In 
this figure, the two parameters are the isotropy of the light dis- 
tribution ySi4 and the fraction of PMT hits within the prompt 
light time window. 
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TABLE X: Summary of backgrounds. The second column gives the 
event rate in the SNO detector D 2 volume (r=600.5 cm) when an 
average production rate was used. The last column gives the esti- 
mated number of events remaining in the 391-day salt phase data 
set after cuts. The internal neutron and y-ray backgrounds are deter- 
mined from independent information and constrained in the analysis. 
The external-source neutrons are extracted along with the signal esti- 
mated from the energy-unconstrained signal extraction fit. For back- 
grounds for which only an upper limit can be determined, the 68% 
CL upper limit is used as a 1 -standard-deviation uncertainty in the 
error propagation. 
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Counts in 
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Neutrons generated inside D 


2 0: 








2 H photodisintegration [U 


, Th] 






qi t+30.4 
=,J ---'-31.5 


2 H photodisintegration [ 24 


Na] 






10.2 + 2.5 


n from fission [U] 


0.43 




U v" 1 


± 


2 H(a,an) 1 H [Th] 


1.9 


n/<g" 


Thy- 


0.93 ± 0.50 


2 H(a,an) 1 H [ 222 Rn] 


0.80 


n/<g" 


Uy-' 


2.89 ± 0.47 


17 - 18 O(a,n) 20 - 21 Ne [Th] 


0.09 


n/Jg~ 


Thy- 


0.03 ± 0.02 


17 - 18 O(a,n) 20 - 21 Ne [ 222 Rn] 


0.20 


n/Jg~ 


Uy" 1 


0.72 + 0.12 


n from atmospheric v 








15.8! 21 6 3 


24 Na from muons 


0.33 


n y _1 




0.14 ±0.14 


muons in SNO 


11240 


n y -1 




< 1 


muons in rock 


0.14 


n y _1 




0.08 ± 0.01 


v~ e "ccp" 


0.03 


n y _1 




0.01 ± 0.01 


~v~ e "ccd" 


1.43 


n y _1 




0.6 ±0.1 


v~ e "ncd"-reactor 


3.24 


n y _1 




1.4 ±0.3 


V e "ncd"-terrestrial 


1.2 


n y -1 




0.5 ±0.1 


CNOv 


1.0 


n y -1 




0.4 ± 0.4 


Total internal-source neutrons 






125 l +iLi 


y-ray s generated uniformly inside D 2 0: 






y from fission [U] 


0.04 


7A<g 


Uf 1 
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y from atmospheric v 








o 9+4.6 
J * Z -4.4 


Total internal-source y-rays 








q 9+4.6 
:, - Z -4 4 



Decays of spallation products throughout D 2 0: 
16 N following muons 16 Ny~' 
Other spallation 1.2 A Zy~ [ 



< 1.3 

< 0.8 



Cherenkov events from radioactivity inside D 2 0: 

ySy decays (U,Th, 24 Na) 3.6^° 

Backgrounds produced outside D 2 0: 

Externally generated neutrons (from fit) 128.5±42.4 

/3y decays (U, Th) in AV, H 2 0, PMTs < 1 8.5 

Instrumental contamination <3 

Isotropic acrylic vessel events < 6.55 



One key point regarding the bifurcated analysis is that both 
sets of cuts must be sensitive to each background class in or- 
der for that class to be included in the background estimate. 
This is true of all known instrumental backgrounds, except for 
the isotropic events from the acrylic vessel. For this reason a 
separate estimate of this background is required. 
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FIG. 16: The distribution of high-level cut parameters for instrumen- 
tal backgrounds and neutrino candidates. The two cut parameters are 
the isotropy /J14 and the fraction of PMT hits within the prompt time 
window. The 16 N calibration source was used to generate the sample 
"neutrino candidate" events and thereby establish the signal window 
for Cherenkov events and to calibrate the cut efficiency. 



B. Isotropic Acrylic Vessel Background 

A class of background events of uncertain origin was iden- 
tified in the D2O phase. These events reconstruct near the 
acrylic vessel and are characterized by a nearly isotropic 
light distribution. They may result from triboluminescence 
as stresses are relieved in acrylic. In the analysis of the 
data from the D2O phase, the events were removed by the 
isotropy and fiducial volume cuts. However, the isotropy cut 
has been relaxed for the salt phase because it would elimi- 
nate too many neutrons, whose Cherenkov light distribution 
is more isotropic in the presence of salt. 

Two independent analyses were performed to estimate the 
background remaining after the fiducial volume selection cut 
of R <550 cm was applied. In the first analysis the isotropy 
distribution of the salt data was fit to the expected shapes for 
neutrons and electrons. The difference between this fit result 
and the actual signal in the high-isotropy region of the spec- 
trum was attributed to isotropic background events in the data 
set. 

In the second method, a cut was applied on the isotropy 
parameter to the pure D2O data set, where the neutrino signals 
are clearly defined. The difference between the actual number 
of events removed by the cut and the predicted loss of neutrino 
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events yields an estimate of the background in the D2O data 
set. This estimate was then scaled to obtain a limit for the salt 
data. 

The combination of the analyses gives a 68% upper limit of 
6.55 events inside the fiducial volume of R <550 cm and with 
a kinetic energy above 5.5 MeV in the salt data set. The same 
analyses were repeated separately on the day and night data 
sets, yielding a day-night asymmetry, defined as the difference 
between the night and the day signal rates normalized by their 
average, of 0.68±0.31. 



C. Photodisintegration Backgrounds from Internal 
Radioactivity 

The deuteron can be dissociated by a y-ray above the bind- 
ing energy of 2.22 MeV. The j3 - y decays of 208 T1 and 214 Bi 
from the 232 Th and 238 U chains emit y-rays that are above this 
binding energy. Neutrons produced through photodisintegra- 
tion are indistinguishable from those produced by the NC re- 
action, therefore measurement of the levels of backgrounds 
inside the detector is crucial for an accurate measurement of 
the total 8 B flux. A concentration of 3.8 x 10" 15 g Th/gD 2 or 
30 x 10~ 15 gU/gD20 in the heavy water would each contribute 
one photodisintegration neutron per day. 

Two independent approaches were developed to measure 
these backgrounds, which can be classified as ex-situ and in- 
situ techniques. 

Three ex-situ techniques have been developed to assay pre- 
cursor radioisotopes of 208 T1 and 214 Bi in the D2O and the 
H2O. The decays of the parent Rn and Ra isotopes are counted 
in a system external to the SNO detector. Two of these tech- 
niques extract 224 Ra and 226 Ra using beads coated with man- 
ganese oxide (MnO x ) [33] or membranes loaded with hy- 
drous titanium oxide (HTiO) [34]. Radioassays typically in- 
volve the circulation of 500 tonnes of D 2 of which approx- 
imately 400 tonnes is flowed through these media. In the 
MnOx technique, Rn daughters from the Ra decays are iden- 
tified by a spectroscopy. In the HTiO technique, the extracted 
Ra atoms are concentrated and identified by f3 - a coinci- 
dences of their decay products. Because the ingress of long- 
lived (Ti/ 2 = 3.8 d) 222 Rn (e.g., by emanation from materi- 
als or ingress from laboratory air) breaks the equilibrium with 
226 Ra, this background in the water is obtained by degassing 
and cryogenically concentrating the dissolved gas from assays 
of approximately 5 tonnes of water [35]. The 222 Rn decays are 
subsequently counted in a ZnS(Ag) scintillation cell. 

Over the entire salt phase, 16 MnOx and 6 HTiO assays 
were conducted at regular intervals. The addition of salt to the 
heavy water affected the radium assay techniques in two ways: 
1) a somewhat reduced radium extraction efficiency was ob- 
served on MnOx (81% compared to 95% in pure D 2 0) and 2) 
dissolved manganese steadily increased from a constant 2 ppb 
in pure D 2 to 12 ppb at the end of the salt phase, causing 
interference with the measurement of the concentration of ra- 
dium in the HTiO assays. During the salt phase, the man- 
ganese could not be removed by the reverse osmosis method, 
as that would have also removed the salt. Since the MnOx 



TABLE XI: 224 Ra assay results for elements of the D 2 systems, 
compared to the total activity measured in the complete system. 
These measurements do not include the small piping section in the 
detector. 



Element of heavy water system 


Th bkgd. ( 224 Ra/day) 


Ultrafiltration unit (assay system) 




Heat exchanger 




Ultrafiltration unit (main stream) 




Main recirculation pump 


7 9+33 


Process degasser 


75 +31 


Filtration unit (main stream) 


< 36 


Piping 


330 ± 184 


Complete system 


816 ± 165 



technique did not suffer from interference with dissolved man- 
ganese in the D2O, its results were used as the ex-situ mea- 
surement of the 224 Ra concentration. 

A small amount of activity was observed in all the elements 
outside the acrylic vessel and corresponded to a total produc- 
tion rate of 8 16 ± 165 224 Ra atoms per day, with approximately 
a quarter of these produced before the columns. This activity 
caused a correction of 0.65 x 10~ 15 gTh/gD20 to the ex-situ 
measurements. An assay showed no significant change in ac- 
tivity when the salt was added to the D2O and after the salt 
was removed most of the activity in the elements outside the 
acrylic vessel went away. The increase in the salt phase (from 
an upper limit of 16 224 Ra atoms per day in the D2O phase) 
is interpreted as primarily due to a displacement of the equi- 
librium between radium in solution and radium bound to its 
Th sources in the water system, arising from the presence of 
Na + ions in the D2O. The time constant associated with the 
removal of this activity showed that it originated from flowing 
as well as stagnant segments of the water circulation and as- 
say systems. Components of the systems can be isolated and 
the 224 Ra contribution from each component was measured 
by circulating the enclosed water over a MnOx column. Ta- 
ble XI provides a summary of the measured production rate of 
224 Ra in various components in the assay and water circulation 
systems. The amount in the piping outside the AV is the dif- 
ference between the total and the sum of the other individual 
elements. The piping within the AV is substantially shorter 
and thinner than that outside and any contribution from this 
section was assumed negligible. 

During a normal assay, the activity measured on the MnOx 
assay columns is the sum of contributions from the D2O and 
from the piping leading to the columns. The effect of this 
source of Th on the levels in the heavy water target over time 
was modeled, given the source distribution in the water pip- 
ing and the water circulation paths and times. The model di- 
vided the piping external to the acrylic vessel into sections, 
and traced the amount of 224 Ra activity added during the as- 
says. 

The mean level of Th concentration in the D2O target was 
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(1.76+Qgg) x 1(T 15 gTh/gD 2 0, where the quoted uncertainty 
reflects the possible distribution of this background activity. 
Combining this in quadrature with the systematic uncertain- 
ties associated with the MnOx assay technique of ^37% gives 
(1.76 ± 0.44(stat)+° ™(syst)) x 1(T 15 gTh/gD 2 for the ex-situ 
analysis. 

The in-situ technique uses pattern recognition on the 
Cherenkov light distribution to determine the equivalent con- 
centration of 232 Th and 238 U in the water. The decays of 208 T1 
produce a more isotropic light distribution than 214 Bi decays 
because of a more complex decay scheme. In the energy 
window 4 < r eff < 4.5 MeV, decays of 208 T1 and 214 Bi in 
the Th and the U chains are the dominant components of the 
Cherenkov signal. By studying events that reconstruct with 
R <450 cm, 208 T1 and 214 Bi are separated statistically by their 
differences in the light isotropy In this energy interval, 
the Cherenkov light from 214 Bi decays is primarily from the 
direct ground state j3 decay with an endpoint of 3.27 MeV, 
while almost every 208 T1 decay emits a 2.614-MeV 7, accom- 
panied by one or more low-energy y-rays and a f3 with an end- 
point energy up to 1.8 MeV. Therefore 208 T1 decays produce a 
more isotropic light distribution than 214 Bi decays. The in-situ 
technique also has the advantage of providing direct determi- 
nation of the background levels during data taking without any 
assumptions regarding temporal variation between assays. 

The statistical separation of the 208 T1 and 214 Bi decay sig- 
nals using j6i4 is shown in Fig. 17, where the probability den- 
sity functions (PDFs) used in the maximum likelihood anal- 
ysis were determined from MC simulations. Solar neutrino 
and 24 Na signals, which are backgrounds to the in-situ analy- 
sis, were constrained in this analysis (an undistorted 8 B neu- 
trino spectrum was assumed). Some of the background 24 Na 
was produced by neutron activation of the D2O target during 
the deployment of calibration sources. Data runs that were 
known to have significant levels of 24 Na and Rn ingress were 
removed in the run selection process described in Sec. Ill A. 
Additional 24 Na was produced in the D2O in the chimney re- 
gion of the acrylic vessel, in a buffer tank used in the circula- 
tion of D2O, and in the water circulation pipes. These regions 
are not well shielded from fast neutrons and y-rays emitted 
from the rock in the underground laboratory. The residual 
contributions from these sources of 24 Na were tallied, and the 
photodisintegration neutron production rate was found to be 
0.064+0.016 d -1 . The 24 Na background contribution intro- 
duced by water circulation was calculated using a water flow 
model and constrained in the maximum likelihood fit. The /3]4 
PDFs were calibrated by 222 Rn and 24 Na calibration spikes in 
the D2O. The latter has a f3u distribution similar to 208 T1 de- 
cays. 

The amplitudes of the 208 T1 and 214 Bi signals determined 
from the in-situ analysis were converted to the equiva- 
lent concentration of 232 Th and 238 U under secular equi- 
librium. The equivalent concentrations integrated over 
the solar neutrino data set are found to be 0.85^[j 4 2(stat) 
+° : 42 (syst) x lO" 15 gTh/gD 2 and 8.28+°; 83 (stat) +}jj(syst )x 
10~ 15 gU/gD20 respectively. 

Results from the ex-situ and the in-situ analyses of the D2O 
are shown in Fig. 18. They are consistent with each other. The 




FIG. 17: In-situ determination of the low-energy background in the 
DtO. The data points represent low-energy events selected by the 
criteria described in the text. The data isotropy (614) distribution is 
fit to a combination of 208 T1 and 214 Bi distributions. Also shown 
are the 24 Na background and solar neutrino contributions that were 
constrained in the fit. The fit result is shown as the sum histogram. 



ex-situ and the in-situ techniques are independent and their 
systematic uncertainties were independently assessed. To ob- 
tain the best measurement of the equivalent concentration of 
232 Th in the D2O target during the 391 -day live time period, 
the weighted mean of the two techniques is used as the back- 
ground input to the integral neutrino flux measurement de- 
scribed in Sec. XI. The 238 U chain activity is dominated by 
Rn ingress which is highly time dependent. In addition, only 
an upper limit could be obtained for the weighted average of 
the ex-situ measurements because of intermittent backgrounds 
in the radon extraction process. Therefore the in-situ deter- 
mination was used for the estimate of this activity as it pro- 
vides the appropriate weighting by neutrino live time. The 
photodisintegration neutron production rates from the natural 
chains were calculated by MC simulations, which show that 
the equivalent of 3.79 fig 232 Th or 29.85 //g 238 U produces one 
photodisintegration neutron per day in the D2O target. 

The photodisintegration neutron production rates from de- 
cays in the 232 Th and 238 U chains in the current data set are 
summarized in Table XII. To determine the temporal varia- 
tion of the internal backgrounds, the data set was divided into 
four time bins. An in-situ determination was performed in 
each of the time periods, and the equivalent concentrations 
of 232 Th and 238 U were found to be relatively constant. The 
in-situ analysis was repeated for the day and the night data 
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TABLE XII: Photodisintegration neutron production rate from de- 
cays of daughters in the 232 Th and 238 U chains in the D2O. The rates 
from the ex-situ and the in-situ techniques are consistent with each 
other. The total rate from the last row is used for all the solar neutrino 
analyses that do not depend on the solar zenith angle. The day (D) 
and night (N) in-situ results are used as inputs to the day-night solar 
neutrino flux asymmetry analysis. 





2J2 Th 
(n d- 1 ) 


(nd" 1 ) 


in-situ (D) 
in-situ (N) 
in-situ (D+N) 


15 +0 - 23 
22 +016 

-0.16 


32 +0 - 05 

"-^-0.05 

0-26^ 
0.28+°$ 


ex-situ (D+N) 


42 +u - 2j 


< 1 


Total (D+N) 


29+u.is 

"■^-0.18 


0.28!"- 



set, and the results are also summarized in Table XII. In the 
day-night asymmetry measurement of the neutrino flux, the 
internal background asymmetry was determined from the in- 
situ analysis of the day and the night data sets, as this gives 
the proper temporal variation. 

Both the ex-situ and in-situ techniques were also applied 
to the determination of radioactive backgrounds in the H2O. 
Throughout the salt phase, 86 radon, 30 MnOx and 13 HTiO 
radioassays of the H2O were performed. Results from the 
radon assays performed on the same day were averaged. 
The MnOx and HTiO results were consistent with one an- 
other and the weighted average, taking into account the neu- 
trino live time, was used to determine the mean concen- 
tration of radioisotopes. The activities were found to be 
5.2+1 i x 10~ 14 gTh/gH 2 and 20.6^ ° x 10~ 14 gU/gH 2 Q by 



-1.6 



the ex-situ techniques. 

In the in-situ analysis of the H2O background, a monitor- 
ing window for events with 4 < T e ff < 4.5 MeV in the 
H2O region (650 < R < 680 cm) was utilized, and the 
equivalent 232 Th and 238 U concentrations were determined 
by fitting the isotropy distribution. The radioactive back- 
grounds in the H2O deduced from the in-situ technique were 
6.1^ x 10~ 14 gTh/gH 2 and W.l^ 1 / x 10~ 14 gU/gH 2 0, 
which are consistent with the results from the ex-situ tech- 
nique. As is discussed below, the neutron background due to 
radioactivity external to the D2O target is determined in the 
fit of the solar neutrino fluxes (Sec. IX). The photodisintegra- 
tion background arising from activity in the H2O is part of this 
external neutron background. 



D. Cherenkov Backgrounds 

The broad energy resolution of the detector allows a small 
fraction of the /?— y decays in the natural radioactive chains to 
appear in the neutrino data sample, even though their Q values 
are lower than the = 5.5-MeV neutrino analysis thresh- 
old. The number of these internal Cherenkov events originat- 
ing within the D2O target is kept small primarily by ensur- 
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FIG. 18: (a) Thorium and (b) uranium backgrounds (equivalent equi- 
librium concentrations) in the D2O deduced by in-situ and ex-situ 
techniques. The MnO x and HTiO radiochemical assay results, the 
Rn assay results, and the in-situ Cherenkov signal determination of 
the backgrounds are presented for the period of this analysis on the 
left-hand side of frames (a) and (b). The right-hand side shows 
time-integrated averages including an additional sampling system- 
atic uncertainty for the ex-situ measurement. The large 222 Rn excess 
near day 580 is the decay of Rn that was added for calibration pur- 
poses, and was excluded in calculating the mean ex-situ results. The 
weighted average of the ex-situ 222 Rn measurements appear as an up- 
per limit only. This is due to intermittent background appearing in 
the radon extraction process. 



ing low radioactivity levels. Outside the heavy water volume, 
however, the acrylic vessel, the light water, and in particu- 
lar the PMT array and support structure have higher levels 
of radioactivity. Most of these external Cherenkov events are 
eliminated by imposing a 550-cm fiducial volume cut. These 
events can "leak" into the fiducial volume in two ways: y- 
rays can travel unscattered from their external origin inward 
and events whose origin is outside the volume can have a mis- 
reconstructed vertex located inside. 

The internal Cherenkov background was determined from 
MC simulations and calibration with a controlled injection of 
81+4 Bq of 222 Rn (a 'radon spike') into the D 2 target. In the 
analyses of the internal Cherenkov background, the ratio be- 
tween the number of internal Cherenkov events and the num- 
ber of detected photodisintegration neutrons above the neu- 
trino analysis threshold was determined for backgrounds from 
the 232 Th and the 238 U chains. This ratio was then normal- 
ized by the measured number of photodisintegration neutrons 
produced and the neutron detection efficiency in the fiducial 



21 



> 




-5 \ -r 

10 r ^ 

3 4 5 6 7 8 9 10 11 12 13 

7 ff (MeV) 

FIG. 19: Fitted MC prediction to the radon spike data, with the 
steeply falling background Cherenkov spectrum and neutron peaks 
shown separately. 



volume (described above). For 208 T1 decays in the Th chain 
and 24 Na, MC simulations of their Cherenkov signals were 
used. The systematic uncertainties were determined by per- 
forming 10 000 simulated experiments, with the scale and res- 
olution of the energy response, vertex reconstruction, and fi\\ 
drawn from distributions estimated by the respective analy- 
ses. All analysis cuts applied to the neutrino data sample were 
also applied here, and their individual uncertainties included 
in the measurement. For the 2I4 Bi decays in the U chain, the 
energy-differential uncertainties of the Rn spike energy spec- 
trum were first determined by 1000 simulated experiments. 
Each of the simulated experiments assumed a different spatial 
distribution of radioactivity within the detector, constrained 
by the reconstructed position of low-energy events in the neu- 
trino data set. Uncertainties associated with the time variation 
of the detector response over the course of the neutrino data 
set were taken from the calibration analyses. The Cherenkov- 
to-photodisintegration neutron ratio was determined by fitting 
the energy distribution of the Rn spike (with its uncertainties 
determined from the simulated experiments) to the simulated 
Cherenkov background and neutron spectra. The parameters 
which were allowed to float in the fit included the energy 
scale, the energy resolution, and the amplitudes of both the 
Cherenkov events and the associated photodisintegration neu- 
trons. Figure 19 shows the resultant fit. The energy scale and 
resolution uncertainties were both consistent with the analy- 
ses of 16 N and 252 Cf data. The systematic uncertainties on the 
ratio were determined by varying all the parameters over their 
allowed uncertainties, including all covariances. 

At r e ff > 5.5 MeV, the ratio of the number of Cherenkov 



events to detected photodisintegration neutrons was found to 
be 0.011;*$* for 2()8 T1, 0.090! ,; , 24 for 24 Na and 0.053!°;°;' 
for 2I4 Bi, which correspond to a total of 3.6^° observed 
events in the 391 -day data set. Since the multiplicative fac- 
tors used to scale the photodisintegration neutron amplitude 
to the number of internal Cherenkov events are independent 
of the solar zenith angle, the internal Cherenkov backgrounds 
in the day and the night data set are obtained by scaling the 
in-situ results in Table XII. 

The analysis methods for the internal Cherenkov back- 
grounds are not directly applicable to the external Cherenkov 
backgrounds. Outside the fiducial volume, the detector is not 
as well calibrated because of greater optical complexity and 
limited accessibility for the primary optical and energy cal- 
ibration sources. To model the radioactivity in this region, 
Th and U sources were deployed at various locations external 
to the fiducial volume and used to generate volume-weighted 
radial distributions for low-energy backgrounds originating 
from the acrylic vessel, the H2O, and the PMT support struc- 
ture. The radial distributions were utilized since they are rel- 
atively insensitive to Th and U differences. For the acrylic 
vessel and the PMT support structure, the p distributions 
were based on a Th source encapsulated in acrylic. For the 
H2O region, a calibrated 222 Rn spike was used. Events with 
r eff > 4.5 MeV and 1.1 < p < 2.5 in the neutrino data 
set were fit by a maximum likelihood technique with radial 
PDFs generated with the same energy threshold as the source 
data. Figure 20 shows the results of this fit, where the width 
of the band in the figure represents the systematic uncertain- 
ties. These uncertainties include spatial and temporal depen- 
dence of the reconstruction, subtraction of contributions from 
photodisintegration neutrons near the acrylic vessel, and the 
difference between the p distributions of the decays of U and 
Th daughters. Due to limited statistics in the source PDFs, 
the fit to the p distribution for the neutrino data set could 
not be done at the solar neutrino analysis energy threshold of 
r e (j > 5.5 MeV. Therefore, the p distributions for the source 
data with T e s > 5.5 MeV were normalized by the amplitudes 
obtained in the fit to the data with T e g > 4.5 MeV and extrap- 
olated inside the fiducial volume to estimate the number of 
background events in the signal window. This analysis was re- 
peated for the day and night data sets. At the neutrino analysis 
energy threshold, the extrapolated contributions of the back- 
grounds from different detector regions are consistent with 
zero. For the 391-day salt data set, the 68% upper limit for 
the acrylic vessel, H2O, and the PMT background contribu- 
tions are 7, 3, and 1 1 events respectively with a combined up- 
per limit of the total external Cherenkov background of 18.5 
events. The day-night asymmetry of this external background 
is-0.10±0.16. 



E. Other Backgrounds 

In addition to the main contributions stemming from 
deuteron photodisintegration and low-energy Cherenkov 
events from (3 - y decays, there are other sub-dominant back- 
grounds that must be assessed. These backgrounds include 
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FIG. 20: Maximum likelihood fit of the data p distribution (data 
points) to the PDFs constructed from Th or U calibration sources. 
Because the external Cherenkov backgrounds are small for the so- 
lar neutrino analysis threshold (T eff > 5.5 MeV) and fiducial volume 
(p < 0.77), this fit is performed for T e g > 4.5 MeV in the region 
1.1 < p < 2.5 in order to enhance the statistics. Contributions of 
these backgrounds are then obtained from extrapolating the fit re- 
sults to the solar neutrino signal window. The band represents the 
systematic uncertainties in this analysis. 



contributions from a reactions on elements in the water target 
and the construction materials, 238 U fission, cosmic ray spal- 
lation, and reactor and atmospheric neutrinos. 

A small source of neutrons can come from a reactions on 
2 H, n O, and ls O. The as produced in the uranium and thorium 
decay chains. Neutrons can also be produced from the sponta- 
neous fission of 238 U, which has a half-life of (8.2±0.1)x 10 15 
years. The concentration of 238 U has been measured using 
ex-situ HTiO radioassays. Results from the assays indicate a 
negligible contribution of neutrons and y-rays from the spon- 
taneous fission of 238 U in the D2O. 

Backgrounds from atmospheric neutrino interactions were 
estimated with the aid of the NUANCE [36] neutrino Monte 
Carlo simulation package. Atmospheric neutrino interactions 
can contribute to the production of neutrons without other en- 
ergy deposits to tag the event, or via the production of un- 
tagged photons from the de-excitation of I6 from neutral- 
current neutrino-nucleon scattering. The NUANCE simula- 
tion provides a comprehensive estimate of various neutrino 
interactions and includes final state intranuclear interactions. 
The estimated contribution to this data set from such events 
passing all selection cuts is 15.8^|g 3 neutrons and 3.2^-JJ y- 



rays. 

Events that produce two or more neutrons within short 
time intervals can serve as a test of the background contri- 
butions from atmospheric neutrinos, spontaneous fission, and 
2 H(«,2n)'H. Analysis of burst data taken from the full data 
set compared to Monte Carlo predictions show an excess of 
high-multiplicity bursts in the neutrino data. A burst is de- 
fined as two or more events, passing all neutrino selection cri- 
teria, that occur within a 50 ms time interval. A likelihood 
calculation indicates that, if the MC accurately describes the 
data, the probability of obtaining a worse likelihood is 1.6%. 
Uncertainties in hadron transport and intranuclear reactions 
for atmospheric neutrino interactions are considered the likely 
cause of the deficit of high-multiplicity events in the MC pre- 
diction. In estimating the uncertainty associated with the sin- 
gle neutron events from atmospheric neutrino interactions, the 
upper uncertainty has been conservatively taken to encompass 
the difference between the data and the Monte Carlo predic- 
tion. 

Neutrons and y-rays produced at the acrylic vessel and in 
the light water can propagate into the fiducial volume. Dur- 
ing construction of the acrylic vessel, Rn progeny accumu- 
lated on its surfaces. These daughters can initiate (a,n) reac- 
tions on 13 C, n O, and 18 0. External y-rays originating from 
(a, ny) and (a, py) processes and radioactivity in the construc- 
tion material of the detector and the rock cavity can enter the 
D2O target and photodisintegrate deuterons. The radial dis- 
tribution of these neutron sources differs from those of the 
NC signal and of photodisintegration neutrons produced from 
radioactivity in the D2O target. The enhanced neutron detec- 
tion efficiency of the salt phase makes it possible to extract 
the external-source neutron contribution in the neutrino signal 
window by including an additional radial distribution function 
in the statistical analysis of the solar neutrino flux. Details of 
the extraction of this background can be found in Sec. IX. Ad- 
ditional tests, including direct counting of the a activity on the 
surface of the acrylic and the search for coincident events gen- 
erated by specific nuclear reactions associated with the (a,n) 
reaction, were performed. The main source of coincidence 
events is the e + e~ pairs from the excited state at 6.05-MeV in 
16 0, in coincidence with a fraction of the neutrons produced 
by n C(a,n) 16 0. A weaker source is the two neutrons from 
2 H(n, 2n)'H induced by fast neutrons from (a,n) reactions. 

The results of these measurements are limited by statistics. 
The sum of this a-induced neutron background and the photo- 
disintegration neutron background produced by radioactivity 
in the H2O and the AV is consistent with the results from the 
radial fit technique. 

At the depth of SNO, only neutrinos and muons from cos- 
mic rays survive. Neutrons and other cosmogenic activity are 
produced from muon capture, muon electrodisintegration, and 
from muon nuclear spallation. The 20-second veto following 
a tagged muon event removes most subsequent activity. How- 
ever, longer-lived spallation products, such as 16 N, can survive 
this cut. A limit of less than 1.3 16 N events (68% CL) was de- 
termined by analysis of long-lived activity present within 50 
seconds after a muon event. Activity from other spallation 
products is estimated to contribute less than one event in the 
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data set. 

A small number of neutrons can also be created by v e reac- 
tions from nuclear reactors. The estimate of this background 
is based on the average power output of all commercial reac- 
tors within 500 km of SNO and an average reactor v e spec- 
trum. Electron antineutrinos from terrestrial radioisotopes do 
not contribute significantly to the background because their 
energies are below threshold for the CC reaction on 2 H ("ccd") 
and can only initiate CC reactions on 'H ("ccp") and NC re- 
actions on 2 H ("ncd"). The total estimated background contri- 
butions from reactor and terrestrial v e interactions are 2.0±0.3 
and 0.5 ±0.1 detected neutrons, respectively. 

SNO is slightly sensitive to solar CNO neutrinos generated 
by the electron capture decay of 15 and 17 F, and this con- 
tributes 0.4 ± 0.4 neutrons [37]. 

None of the backgrounds discussed in this section depend 
on the solar zenith angle. Thus, the contribution of these back- 
grounds to the day and the night data sets can be determined 
by normalizing to the respective live times. 



VIII. SYSTEMATICS ASSOCIATED WITH THE 
DAY-NIGHT ASYMMETRY MEASUREMENT 

Differences in day and night neutrino fluxes are a predic- 
tion of matter-enhanced neutrino oscillations. Day-night re- 
sults are reported as asymmetry ratios in the measured neu- 
trino fluxes. The asymmetry ratio for a flux is defined as 
A = 2(<f>pf - <Pd)I{<Pn + 4>d)- An advantage of the asymme- 
try ratio is that most systematics cancel and only systematic 
effects that scale day and night fluxes by different factors need 
to be considered. 

Day-night systematics can be divided into four general 
classes. Diurnal systematics are variations in detector re- 
sponse over a 24-hour timescale, such as might be caused by 
diurnal changes in the laboratory environment. Because the 
SNO detector is located far underground, it is isolated from 
many diurnal effects. Day-night differences in detector re- 
sponse are therefore not expected, but limits must be placed 
on their size. Directional systematics arise because the SNO 
detector is not completely spherically symmetric, and because 
the directions of electrons from CC or ES neutrino interactions 
are correlated with the time of day. ES events in particular 
are highly directional, and so ES events at night will prefer- 
entially illuminate the upper half of the detector, while during 
the day they illuminate the bottom half. If there are differences 
in the up-down response of the detector, these directional dif- 
ferences can create effective differences in the day and night 
rates. Directional systematics are expected to be important for 
ES events, but are greatly suppressed for CC events, which 
have only a weak directional correlation with the direction 
from the Sun. Directional systematics do not produce day- 
night systematic effects for directionally isotropic events such 
as neutrons or backgrounds. Miscellaneous systematics in- 
clude possible day-night differences in cut acceptance, uncer- 
tainties in the live time calculation, and long-term variations in 
detector response. These can produce differences in the time- 
averaged day and night detector responses if the long-term 



TABLE XIII: Day-night differences for selected muon-induced neu- 
tron distributions. 





Asymmetry 


Quantity 


Night value Day value ratio (%) 



Event rate (day" 1 ) 3.50 ± 0.09 3.56 ±0.10 -1.83 ±3.81 
Mean energy (Me V) 5.78 ± 0.06 5.67 ± 0.06 1.96 ± 1.49 
Energy width (MeV) 1.47 ± 0.06 1.47 ± 0.06 0.00 ±4.69 
Mean isotropy (J3 lA ) 0.31 1 ± 0.004 0.312 ± 0.005 -0.24 ± 2.03 
Capture time (ms) 4.9 ± 0.3 4.6 ± 0.3 7.27 ± 8.36 



variations are correlated with the seasonal variations in the 
day-night live time exposure. Background systematics reflect 
uncertainties in the magnitude and day-night asymmetries of 
the various background sources. The following sections de- 
scribe how each class of systematics is measured. 



A. Diurnal Systematics 

Diurnal systematics are variations in detector response with 
a 24-hour periodicity. A strictly hypothetical example would 
be diurnal variations in the laboratory's temperature, which 
in principle could change the response of the SNO electron- 
ics if there were uncompensated temperature dependencies. 
Because the timescale for diurnal variations is much shorter 
than the average interval between calibrations, diurnal stabil- 
ity must be evaluated using classes of events that are con- 
tinually present in the detector. These primarily consist of 
secondary neutrons produced in the D2O by through-going 
muons and Cherenkov events from low-energy B-y decays. 

Muons traversing the acrylic vessel can produce secondary 
neutrons by a variety of processes, including photodisintegra- 
tion of deuterons and nuclear spallation. These neutrons are 
produced uniformly throughout the D2O and at an essentially 
constant rate on a diurnal timescale, and so mimic the solar 
NC signal in spatial and temporal distributions. Small diur- 
nal and seasonal variations of a few percent in the predicted 
muon rate, expected as a result of variations in the scale height 
of the atmosphere, are neglected. Secondary neutrons from 
muons are selected by identifying bursts of events inside the 
D2O that occur within a time window of 50 //s - 20 ms follow- 
ing a tagged muon event. Spallation events can occasionally 
produce very large bursts of neutrons. To prevent such bursts 
from biasing the spatial or temporal distributions of the neu- 
trons, a multiplicity selection requiring less than 15 events in 
the bursts is applied to the neutron selection. 

Table XIII shows the mean day and night values for various 
distributions of muon-induced secondary neutrons. The event 
rate, mean event energy, width of the neutrons' energy distri- 
bution, mean value of the Bu isotropy parameter, and neutron 
capture time are all consistent between day and night. The 
left side of Fig. 21 shows the day-night asymmetry on each 
measured quantity. 

While muon-induced secondary neutrons demonstrate the 
diurnal stability of the detector, they are statistically limited 
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TABLE XIV: Combined day-night asymmetries for various back- 
ground regions. The asymmetry ratio is formed by combining the 
measured day-night asymmetries for each run with both day and 
night data. 



Background region 


Combined asymmetry ratio 


D 2 


-0.045 ± 0.070 


AV hot spot 


-0.022 ± 0.048 


Shell around AV 


-0.023 ± 0.022 


H 2 


0.003 ± 0.016 


PMTs 


0.009 ±0.011 



by the low muon rate (~ 3/hour). Better limits on most pa- 
rameters may be obtained from studies of intrinsic detector 
radioactivity. These include Cherenkov events from /5-y de- 
cays of low-energy radioactivity in the D2O, acrylic, H2O, or 
PMTs. 

A localized region of higher-than-average background ra- 
dioactivity was discovered on the upper surface of the acrylic 
vessel. The origin of this "AV hot spot" is uncertain, but it was 
presumably introduced by radioactive contamination during 
construction. If composed of Th, the total amount is approxi- 
mately 10 ;Ug, and is not a significant neutron source given the 
location. However, this hot spot provides an excellent check 
of position reconstruction for a point source of events. Com- 
parisons of the reconstructed day and night positions of the 
AV hot spot show that its position is stable to within ±0.3% 
in radius between night and day, and indicate that the vertex 
resolution for a point source differs by no more than ±1 .26 cm 
between night and day. 

Because low-energy background events have steeply falling 
energy spectra near and below the analysis threshold, small 
variations in energy scale or energy resolution produce large 
changes in the observed rate inside a low-energy window. 
Measurements of diurnal rate stability for low-energy radioac- 
tivity can thus be used to limit diurnal variations in energy 
scale and energy resolution. 

In order for this procedure to work, the actual level of ra- 
dioactivity must be diurnally constant. Radioactivity inside 
the PMTs or acrylic is immobile, and presumed to be constant 
except for possible slow decay. Radioactivity in the D2O and 
H2O can fluctuate over time due to radon ingress, as is seen in 
the in-situ radioactivity measurements, but the timescale for 
these changes is generally long compared to the 12-hour dif- 
ferences being sought in this analysis. These longer-term vari- 
ations in radioactivity can be mitigated by calculating a day- 
night asymmetry for each individual run that includes both 
day and night live time (implicitly assuming that the radioac- 
tivity is constant over the several hour duration of a typical 
run). Assuming the real radioactivity level does not vary sig- 
nificantly over this short time period, any observed rate varia- 
tion would indicate a diurnal difference in detector response. 
The run-by-run asymmetries can be combined in a weighted 
average using a maximum likelihood technique to determine 
an overall limit on the rate asymmetry of each source. 

The right side of Fig. 21 shows the combined run-by-run 
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FIG. 21: Day-night asymmetries of selected muon-induced neutron 
properties and background rates. The day-night asymmetries for the 
muon-induced neutron properties are calculated for events in the full 
day and night data sets. The day-night asymmetries of the back- 
ground rates are determined by combining the calculated asymme- 
try of each individual data run. Each data run lasted for less than 
24 hours. 



asymmetries for radioactivity in five different regions of the 
detector: in the AV hot spot, in a spherical shell around the 
AV (excluding the hot spot), in the light and heavy water, and 
near the PMTs. Events used to calculate these asymmetries 
were selected from an energy region corresponding to A^, lompt 
between 25 and 40 (approximately 4.1 < T e g < 6.2 MeV). 

Table XIV shows the asymmetries for each of these back- 
ground regions. All are consistent with no asymmetry. The 
analysis was repeated, this time calculating day-night asym- 
metry ratios for sets of 24-hour periods instead of for each 
run, and similar results were obtained. The insensitivity to the 
length of the normalization period of the asymmetry, and the 
fact that asymmetries in all regions of the detector are consis- 
tent with zero, make it very unlikely there are diurnal varia- 
tions in the detector response. 

While the rate asymmetries on radioactivity in the H2O and 
near the PMTs have the smallest uncertainties, these regions 
are far from the fiducial volume of the neutrino analysis. A 
more representative approach is to base the energy scale and 
resolution variation estimates on the rate of radioactivity in the 
spherical shell around the AV. This region has adequate statis- 
tics, is at a similar radius compared to neutrino events, and 
includes a significant fraction of acrylic radioactivity which 
should be diurnally constant. (This region also contains events 
from thin shells of D2O and H2O adjoining the AV.) Allowing 
the rate asymmetry for this region to vary within its measured 
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limits, and assuming that the entire variation is due to changes 
in energy scale, the diurnal energy scale variation is limited 
to ±0.4%. Assuming instead that the entire variation is due 
to changes in energy resolution, the diurnal energy resolution 
variation can be limited to ±0.5%. These results are consis- 
tent with limits derived from the other regions of radioactivity 
and with the more direct measurements from muon-induced 
neutrons given in Table XIII. 

In-situ radioactive background measurements can also 
place limits on diurnal variations in event isotropy. Compar- 
ison of the mean day and night fiu values for events in the 
spherical shell around the AV results in the limit |/?i4,ni g ht - 
j8i4,day| < 0.006. This limit is comparable to the indepen- 
dent limit of 0.0064 derived from muon-induced neutrons (Ta- 
ble XIII), and the two limits may be combined to limit the 
diurnal shift in the mean value to < 0.0043. 



B. Directional Systematics 

As described near the start of this section, variations in de- 
tector response with the direction of the event can produce 
a day-night systematic for neutrino signals that have differ- 
ent directional distributions for night and day. The PDFs for 
the direction of electrons produced by CC and ES interactions 
with respect to the Sun's direction are shown in Fig. 23(b). 
Because the y-rays emitted by neutron captures have random 
directions, detector asymmetries produce no day-night varia- 
tions for NC events. 

Variations in detector response with direction are mea- 
sured with 6.13-MeV y-rays from the 16 N source. Events 
are grouped by bins in cos 6 and (f> in detector coordinates. 
For each bin, the means of each of the eight event variables 
listed in Table XV are calculated. 16 N calibration runs at dif- 
ferent positions in the heavy water are combined in a volume- 
weighted average. The result is a map of the detector response 
versus directional bin for each systematic. 

Monte Carlo calculations predict how many CC and ES 
events fall in each directional bin for SNO's live time expo- 
sure for both night and day. By convolving the MC predic- 
tion for the directional distribution in both time bins with the 
measured detector asymmetry, the mean value of a detector 
response function is calculated for the night and day data sets. 
The volume-weighted difference between the night and day 
values gives a measure of the directional systematic for each 
neutrino signal. 

Table XV gives the limits on the effective day-night differ- 
ence in detector response for CC and ES events resulting from 
directional asymmetries. 



C. Miscellaneous Systematics 

The uncertainties on the day and night live times are calcu- 
lated to be ±0.021%. These result in a 0.03% uncertainty on 
the day-night asymmetry. 

The cuts used to reject instrumental backgrounds are de- 
scribed in Sec. IIIB. The time variability in the fraction of 



TABLE XV: Directional day-night systematics for CC and ES events. 
The day-night effect of directional systematics for NC and back- 
ground events is zero. 



Systematic 



Limit for CC events Limit for ES events 



Energy scale ±0.09% 

Energy resolution ±0.13% 

Radial shift ±0.02% 

Vertex resolution ±0.13% 

Angular resolution, a M ±1.4% 

Angular resolution, p M ±0.7% 

Angular resolution, J3 S ±0.6% 

Isotropy ±0.09% 



±0.79% 
±1.3% 
±0.15% 
±1.4% 
±11.6% 
±6.1% 
±5.2% 
±0.82% 



good events removed can be measured by applying the cuts to 
muon-induced secondary neutrons. The day-night asymmetry 
of the fraction of good events removed by the cuts is measured 
to be A = 0.18 ± 0.33. Note that this is an asymmetry on a 
very small loss fraction. 

In addition to diurnal variations in detector response (see 
Sec. VIII A), variations on longer timescales could exist. Such 
variations can indirectly introduce differences in the day-night 
detector response if they correlate with seasonal variations 
in the day-night live time exposure. For example, if the 
energy scale were slightly mis-calibrated during the winter 
months, this miscalibration would affect the night data set 
more than the day data, since during the winter more night 
data is collected than day data due to the seasonal variation in 
the lengths of night and day. The effects of these variations 
can be constrained by constructing worst-case models that 
systematically over-estimate the measured response in sum- 
mer and under-estimate it during the winter, or vice versa. 
These worst-case models were constructed for energy scale 
and isotropy variations using the regular calibration points 
taken with the 16 N source and their effect is shown in Ta- 
bles XXV and XXVII as long term variations in these pa- 
rameters. 



D. Day-Night Background Systematics 

Uncertainties in backgrounds to the neutrino signals can 
produce day-night systematics in two ways. First, there 
may be uncertainty in the day-night variation of the back- 
ground (i.e., an uncertainty on the rate asymmetry A of the 
background). This uncertainty will differentially affect the 
amounts of background subtracted from the night and day data 
sets. Second, the uncertainty on the total amount of back- 
ground (night+day) results in an uncertainty on the average 
neutrino flux, which enters into the denominator of the day- 
night ratio. Generally, the uncertainty on the asymmetry is 
a larger effect than the uncertainty on the amount of back- 
ground, but both contribute. 

Details of the day-night calculations of the background to- 
tals have been given previously. The reader is referred to the 
discussion of the individual backgrounds in Sec. VII. 
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IX. NEUTRINO SIGNAL DECOMPOSITION 

Extraction of the electron energy spectrum, flux and day- 
night asymmetry of the 8 B solar neutrinos is carried out via 
an extended maximum likelihood fit of the event variables in 
the 4722-event data set by MC generated PDFs. The data vari- 
ables utilized are T e s, p, cos G Q , and p 1 ^, and their distributions 
are shown in Figs. 22 and 23. 

In order to obtain the electron energy spectra of CC and 
ES interactions, PDFs were created for T e ff intervals which 
spanned the range from 5.5 MeV to 13.5 MeV in 0.5 MeV 
steps. For T e ft values between 13.5 and 20 MeV, a single bin 
was used. Minor adjustments were applied to the PDFs to take 
into account signal loss due to instrumental cuts not modeled 
by the simulation (Fig. 2). PDF normalizations for CC and ES 
components were separately allowed to vary in each r e ff bin to 
obtain model-independent spectra. Only the overall normal- 
izations of NC and external neutron components were allowed 
to vary since their T e ft spectra are simply determined by the 
energy release following neutron capture on 35 C1 or 2 H. 

The parameter fiu is correlated with r e ff and, to a lesser ex- 
tent, with p. Similarly, cos 9 Q is weakly correlated with r e ff 
and p. These dependencies were taken into account through 
the use of a multi-dimensional PDF P(T e g,J3u,p, cos Q Q ) fac- 
torized as follows: 

P(T eS ,p u ,p, cos 6 e ) = P(T en ;/3 14 ,p)xP(cos0 Q \T eS ,p), (12) 

where the first factor is just the 3-dimensional PDF for the 
variables 7^, j3\4, and p, while the second factor is the condi- 
tional PDF for cos 9 Q given T e g and p. This approach explic- 
itly preserves all correlations between the four relevant param- 
eters with the exception of the correlation between cos 9 e and 
j8i4, which is assumed to be linked only through T e ff and p. 

To confirm our understanding of these correlations and to 
verify the results with an independent approach, extraction 
of solar neutrino results has also been performed with the 3- 
dimensional PDF in Eq. 12 further reduced to the following 
formulation: 

P(T eS ,p u ,p, cos ) = P(T eS ,p u ) x P(cos e ) x P(p). (13) 

The signal extraction procedure was applied to 100 simulated 
data sets each generated to simulate the expected character- 
istics of the data. It was found that the parameterization in 
Eq. (13), which ignored the correlations with p and cos Q , re- 
sulted in a small bias compared with the approach in Eq. (12). 
As a further cross-check, both approaches were applied to the 
data and yielded results that were consistent with this inter- 
pretation. After applying corrections for the expected bias, the 
results were found to be nearly identical. The average differ- 
ence in the extracted CC signal between the two factorization 
approaches, after corrections, is less than 1%. 

In the flux measurement, a "constrained" fit, in which the 
CC energy spectral shape is fixed as an undistorted 8 B spec- 
trum [38], and an "unconstrained" fit, in which this constraint 
is removed, are carried out. The constrained analysis is useful 
in testing the null hypothesis of neutrino flavor transforma- 
tion under the assumption of an undistorted 8 B solar neutrino 



> 




p 



FIG. 22: (a) T e g and (b) p distributions for CC, ES, NC and external 
neutron events. Where internal and external neutron distributions 
are identical the distribution is simply labeled neutrons. Note that 
the distribution normalizations are arbitrary and chosen to allow the 
shape differences to be seen clearly. The CC energy spectrum shape 
corresponds to an undistorted 8 B model. 



spectrum. The unconstrained analysis has the advantage that 
this model assumption is removed. 

Compared to the pure D2O phase, the addition of salt in- 
creases the sensitivity at large radius to neutron capture (see 
Fig. 14) making it possible to detect background neutrons 
originating at or near the AV. As shown in Fig. 22(b), the dif- 
fering radial profiles allow external neutrons to be separated 
from neutrons generated within the D2O volume. In the pure 
D2O phase analysis [4], the amplitudes of the photodisintegra- 
tion neutron backgrounds from internal radioactivity in the AV 
and the H2O were fixed in neutrino signal decomposition. In 
the current analysis, the p distribution provides discrimination 
between the external neutron and the signal flux contributions. 
The external neutron component is measured simultaneously 
with the flux signals from the likelihood fit. 

In the salt phase analysis, the amplitude of the p PDF of 
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FIG. 23: (a) /?, 4 and (b) cos6» G distributions for CC, ES, NC and 
external neutron events. Where internal and external neutron distri- 
butions are identical the distribution is simply labeled neutrons. Note 
that the distribution normalizations are arbitrary and chosen to allow 
the shape differences to be seen clearly. 



the external source neutrons is allowed to vary in the maxi- 
mum likelihood fit. If the production rates of these external 
neutron backgrounds were constant, they would contribute an 
amount smaller than the systematic uncertainties reported in 
the earlier pure D2O phase analysis. 

Results from the measurements of the radioactive back- 
grounds (Sec. VII) were incorporated in the extraction of so- 
lar neutrino signals. The internal photodisintegration neutron 
backgrounds are subtracted from the fitted NC event total. A 
separate PDF for internal y-rays from atmospheric neutrino 
interactions was included with a normalization fixed to 3.2 
events (in accordance with Table X) and with a shape based 
on events generated by the I6 N calibration source because of 
their similarity. Other backgrounds, listed in Table XVI, were 
treated as systematic uncertainties applied to the appropriate 
T e ff intervals for CC and NC components following signal ex- 
traction. The distinctive dependence of ES events on cos 
was assumed to reduce the effect of these backgrounds on the 



TABLE XVI: Additional background contributions treated as sys- 
tematic uncertainties applied to appropriate T e g intervals in the CC 
and NC spectra. Note that the r cff range for the first bin is 5.5-6.0 
MeV. 



Systematic parameter 
Internal y bkgd. 



ltr limits 



± 2.2 events 



Internal Cherenkov bkgd. - Tl 
Internal Cherenkov bkgd. - Bi 
Internal Cherenkov bkgd. - Na 
External Cherenkov bkgd.- PMT 
External Cherenkov bkgd.- H2O 



<0.53 events (1st CC bin) 

< 2.29 events (1st CC bin) 

< 0.90 events (1st CC bin) 

< 1 1.0 events (1st CC bin) 
<3.0 events (1st CC bin) 



External Cherenkov bkgd. 
AV events 

Instrumental bkgd. 



AV < 7.0 events (1st CC bin) 
< events (1st CC bin) 

< S|5 events (T cfr >6.5 MeV) 

< 3 events (across T c r range) 



ES results to negligible levels. 

Day-night asymmetries are determined for both the con- 
strained and unconstrained fit cases. As described below, 
asymmetries are calculated both allowing all signal fluxes to 
vary and with the NC asymmetry constrained to zero. 

The maximum likelihood fit returns the estimated numbers 
of CC, NC, ES, and external neutron events, with statisti- 
cal uncertainties. Final fluxes are determined by normaliz- 
ing to the solar-model prediction calculated by MC simulation 
with several correction factors, summarized in Table XVII, 
applied. Additional corrections are applied to adjust the pre- 
dicted event rates as a result of CC interactions on n O, ls O, 
23 Na, and 35 ' 37 C1, which are not modeled in the MC simula- 
tion. The MC model uses the Effective Field Theory (EFT) 
calculation in Ref. [39] to obtain the neutrino-deuteron cross 
sections. This EFT calculation is normalized to the standard 
potential-model calculation [40, 41] by fixing the two-body 
axial exchange-current counter term Lj a- Corrections are ap- 
plied to the flux results to account for small differences in the 
choices of the axial coupling constant gA [42] that were used 
in the theoretical calculations, and for normalization to the im- 
proved potential-model calculation [41] (the EFT calculation 
was done with Lj a = 5.6 fm 3 ). The values of the fundamen- 
tal constants [43] used in the present calculations of the cross 
sections, together with correction factors derived to bring the 
cross section calculations in accordance with those constants, 
are listed in Table XVII. Radiative corrections included in 
[41] were taken out, and the corrections of Kurylov et al. [44], 
parametrized as follows, were applied to the CC, NC and ES 
cross sections: 



wee 
wnc 
wes 



1.0318 
1.0154 



7.45 x 10~ 4 £ e 



4.72 x 10~ 6 £? 



= 0.9764-7.81 x 10" 4 r eff 

-1.31 x 10" 4 r e 2 ff + 3.64 x io^ 6 r ( 



(14) 
(15) 

(16) 



where E e is the true total energy of the electron. Table XVII 
also includes the corrections to the CC, ES, and NC fluxes 
from CC interactions on oxygen, sodium, and chlorine, by 
which, together with the radiative corrections, the EFT and 
ES cross sections were multiplied before comparing with the 
data to extract the fluxes. 
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TABLE XVII: Fundamental constants used in the MC and post-fit 
correction factors to the neutrino flux in different reaction channels. 
The radiative correction functions are discussed in the text. 



Constants used 


8a 


1.2670(30) 




G F 


1.16639(1) 10~ 5 GeV- 2 


L\,A 


4.0 fm 3 




sin 2 9 W (MS) 


0.23113(15) 




Correction 


CC ES 


NC 


CC on 0, Na and CI 


1.0081 1.0 


1.0 


Radiative corrections 


<^CC W ES 


^>NC 



Data 
Fit result 
Neutrons 
CC 
ES 

External neutrons 
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FIG. 24: Data spectrum with statistical uncertainties. Included 
are MC spectra for neutron, CC, ES, and external neutron distribu- 
tions. Note that an undistorted 8 B spectral shape has been assumed 
and each MC contribution has been normalized to the number of cor- 
responding fit events measured by the energy-constrained signal ex- 
traction. 



The flux measurements are presented in terms of the CC, 
NC, and ES signals and as the flux of electron type (<p e ) and 
non-electron type (^ T ) active neutrinos. In general, the ef- 
fects of systematic uncertainties are evaluated by re-fitting the 
data after perturbing the model PDFs by the lcr uncertainties 
determined from calibration and background measurements. 
The differences between the nominal flux fit values and those 
obtained with the systematically shifted PDFs are quoted as 
68% C.L. uncertainties. 

Figure 24 shows the energy spectrum, with statistical uncer- 
tainties, of the data that passes all selection cuts. Included in 
the figure are MC generated spectra for CC, ES, internal neu- 
tron, external neutron components, and their sum. Note that 
the MC generated distributions correspond to an undistorted 
8 B neutrino spectrum and each spectrum has been normalized 
to correspond to the total number of fit events for the given 
component as extracted by the energy-constrained fit. 
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FIG. 25: (a) Extracted CC T e ff spectrum. Systematic uncertainties 
have been combined in quadrature and include only the effect of PDF 
shape change, (b) Extracted CC spectrum with statistical error 
bars compared to predictions for an undistorted 8 B shape with com- 
bined systematic uncertainties, including both shape and acceptance 
components. The systematic error bands represent the fraction of the 
total uncertainty attributable to the given quantity. 



X. SPECTRUM 

Figure 25(a) shows the CC energy spectrum extracted from 
the energy-unconstrained fit. The PDF shape change compo- 
nent (see below) of the systematic uncertainties is added in 
quadrature with the statistical error to provide a combined er- 
ror for each bin. The analogous ES spectrum is presented in 
Fig. 26(a). 

Systematic uncertainties on the extracted CC energy spec- 
trum are calculated by separately varying the PDFs according 
to the estimated lcr uncertainties on the detector parameters. 
Signal extraction is then repeated and the differences between 
the nominal fit values and shifted PDF fit values taken as the 
spectrum systematic uncertainties. 

Systematic uncertainties on the spectrum are divided into 
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FIG. 26: (a) Extracted ES T e g spectrum. Systematic uncertainties 
have been combined in quadrature and include only the effect of PDF 
shape change, (b) Extracted ES r c ft spectrum with statistical error 
bars compared to predictions for an undistorted 8 B shape with com- 
bined systematic uncertainties, including both shape and acceptance 
components. The systematic error bands represent the fraction of the 
total uncertainty attributable to the given quantity. 



TABLE XVIII: CC spectrum combined shape and acceptance sys- 
tematic uncertainties for three sample energy ranges. An undistorted 
8 B shape has been assumed. Note that the full table is given in the 
Appendix. 



Source 



CC differential 
(5.5-6.0MeV) (7.0-7. 



uncertainty (%) 

5 MeV) (9.0-9.5 MeV) 



Energy scale (const.) 


-8.0, 


9.5 


0.6, 


0.0 


4.7, 


-5.0 


Energy scale (E dep.) 


-1.5, 


1.6 


-0.4, 


0.6 


1.3, 


-0.4 


Energy radial bias 


-5.7, 


6.1 


0.1. 


-0.1 


3.8, 


-2.9 


Energy resolution 


4.9, 


-4.9 


-0.1, 


0.1 


-2.6, 


2.6 


1 4 scale 


7.1, 


-8.4 


4.4. 


-4.6 


1.1, 


-1.7 


/?14 width 


-0.9, 


-0.1 


0.0, 


-0.2 


0.3, 


-0.2 


Radial scale (const.) 


-2.6, 


2.5 


-2.4, 


2.6 


-2.5, 


2.8 


Radial scale (E dep.) 


0.2. 


-0.2 


-1.2, 


1.2 


-2.1, 


2.1 


Vertex x 


-0.3. 


0.0 


-0.5, 


0.0 


0.1. 


-0.1 


Vertex y 


-0.2, 


-0.3 


-0.2, 


().() 


0.1, 


0.4 


Vertex z 


-0.4, 


-0.2 


0.2, 


-0.6 


-0.3, 


-0.1 


Vertex resolution 


-0.4, 


0.4 


-0.3, 


0.3 


0.4, 


-0.4 


Angular resolution 


-0.2, 


0.2 


0.1, 


-0.1 


-0.6, 


0.6 


Internal y 


0.4. 


-0.6 


0.1, 


-0.1 


0.0, 


0.0 


Selection efficiency 


-0.2, 


0.2 


-0.1, 


0.2 


-0.1, 


0.2 


Backgrounds 


-8.5, 


0.0 


-0.1, 


().() 


-0.1. 


+0.0 



has been assumed for generation of these uncertainties. The 
complete table of uncertainties for all energy bins is included 
in the Appendix. 

Figure 27 shows the CC spectrum PDF shape systematic 
effects and combined systematics for each of the four energy- 
related systematic uncertainties versus T er y under the assump- 
tion of an undistorted 8 B shape. Figures 25(b) and 26(b) show 
the extracted CC and ES spectra as a functions of T e g, with 
statistical error bars on the data, and combined systematic un- 
certainty error bands centered on the prediction for an undis- 
torted 8 B shape (normalized to the data). Note that the error 
bands correspond to the fraction of the total error, calculated 
as the ratio of the square of the given uncertainty to the square 
of the total, attributed to the indicated source (e.g., f3\4 sys- 
tematic errors). 



XI. SOLAR NEUTRINO FLUX RESULTS 



those caused by changes to the PDF shapes, and those result- 
ing from uncertainties in overall acceptance. The former lead 
to errors in the fitted number of charged current events in each 
spectral bin, while the latter lead to errors in the translation 
of differential event counts into differential neutrino fluxes. 
For example, shifting the energy scale alters both the PDF 
shapes and the acceptance in each bin, while a change in the 
radial scale (and hence the effective fiducial volume) primar- 
ily affects the acceptance. For a particular neutrino spectrum 
model, the fitting uncertainties from the PDFs shapes and the 
uncertainties on the acceptance have correlations with each 
other which must be taken into account when calculating bin- 
by-bin neutrino fluxes. A list of these uncertainties is given 
in Table XVIII for the r eff intervals 5.5 - 6.0 MeV, 7.0 - 7.5 
MeV and 9.0 - 9.5 MeV. Note that an undistorted 8 B shape 



As described in a previous section, a statistical analysis to 
separate the neutrino candidate events into CC, NC, ES, and 
external-source neutrons was performed using an extended 
maximum likelihood technique using the [iu, cos 6> G and p 
distributions. The extended maximum likelihood analysis 
yielded 2176+78 CC, 279±26 ES, 2010+85 NC, and 128+42 
external-source neutron events (recall that 125.4 internal neu- 
trons and 3.5 internal y-ray events have been subtracted off). 
Note that a fixed hep contribution of 9.3 x 10 3 crrT 2 s _I has 
been assumed. The systematic uncertainties on the derived 
fluxes for this "energy-unconstrained" fit are shown in Ta- 
ble XIX. The ySi4, cos f? , and p distributions for the selected 
events are presented in Fig. 28 with only statistical uncertain- 
ties shown. 

For the unconstrained analysis, the quoted CC and ES 
fluxes are the equivalent fluxes of 8 B electron-neutrinos, as- 
suming an undistorted 8 B energy spectral shape, that would 
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FIG. 27: (a) PDF shape change contributions and (b) total contribu- 
tions from energy-related systematic energy uncertainties to the ex- 
tracted CC spectrum assuming an underlying, undistorted 8 B shape. 



produce the same CC and ES event rates above the analysis 
threshold of T e g - 5.5 MeV. For the NC case, the quoted 
flux is the flux of all active neutrino types that would pro- 
duce the same NC rate above the reaction threshold of 2.2 
MeV. The fitted numbers of events give the equivalent 8 B 
fluxes [46, 47] (in units of 10 6 cnrV 1 ) 
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and the ratios of the CC flux to NC and ES respectively are 
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FIG. 28: Distribution of (a) /3i4, (b) cosf? and ( c ) volume-weighted 
radius p. Points with error bars represent data while the MC predic- 
tions for CC, ES, NC + internal and external-source neutron events, 
all scaled to the energy-unconstrained fit results, are as indicated in 
the legend. The dark solid lines represent the summed components. 
The (a) and (b) distributions are for events with T c ff > 5.5 MeV and 
Rflt ^ 550 cm, and are averaged assuming an undistorted 8 B spec- 
trum. The same energy cut has been applied for (c) but events are 
shown out to p < 1.6, where p = 1.0 is the edge of the heavy wa- 
ter volume. The dashed vertical line represents the 550 cm fiducial 
volume cut. 
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TABLE XIX: Systematic uncertainties (%) on fluxes for the energy- 
unconstrained analysis of the salt data set. Note that "const." de- 
notes an energy-independent systematic component and "E dep" an 
energy-dependent part. 



Source 


NC uncert. (%) 


CC uncert. (%) 


ES uncert. (%) 


Energy scale (const.) 


-3.3, 


+3.8 


-0.9, 


+ 1.0 


-1.6, 


+ 1.9 


Energy scale (E dep.) 


-0.1, 


+0.1 


-0.1, 


+0.1 


-0.1, 


+0.1 


Energy radial bias 


-2.0. 


+2.1 


-0.6, 


+0.7 


-1.1, 


+ 1.2 


Energy resolution 


-0.8. 


+0.8 


-0.2. 


+0.2 


-0.7, 


+0.7 


P14 mean (const.) 


-3.6, 


+4.5 


-4.0, 


+3.7 


-1.2, 


+ 1.3 


jS]4 mean (E dep.) 


-0. 1 . 


+0.2 


-0.2, 


+0.0 


-0.0, 


+0. 1 


P14 width 


-0 


+U.U 


-0 2 


+0 2 


-0 2 


+0 2 


Radial scale (const.) 


-3.0. 


+3.3 


-2.6^ 


+2.5 


-2.6! 


+3.0 


Radial scale (E dep.) 


-0.6, 


+0.5 


-0.9, 


+0.8 


-0.7, 


+0.8 


Vertex A' 


-0.0, 


+0.0 


-0.0. 


+0.0 


-0.1, 


+0.1 


Vertex v 


-0.1. 


+0.0 


-0.0, 


+0.0 


-0.1, 


+0.1 


Vertex z 


-0.2, 


+0.2 


-0.1, 


+0.1 


-0.0, 


+0.0 


Vertex resolution 


-0.1. 


+0.1 


-0.1. 


+0.1 


-0.1, 


+0.1 


Angular resolution 


-0.2. 


+0.2 


-0.4, 


+0.4 


-5.1, 


+5.1 


Internal neutron bkgd. 


-1.9, 


+ 1.6 


-0.0, 


+0.0 


-0.0, 


+0.0 


Internal y bkgd. 


-0.1. 


+0.1 


-0.1. 


+0.1 


-0.0, 


+0.0 


Internal Cherenkov bkgd. 


-0.9. 


+0.0 


-0.9. 


+0.0 


-0.0, 


+0.0 


External Cherenkov bkgd. 


-0.2. 


+0.0 


-0.2, 


+0.0 


-0.0, 


+0.0 


Instrumental bkgd. 


-0.4. 


+0.0 


-0.3, 


+0.0 


-0.0, 


+0.0 


Neutron capture eff. 


-2.3, 


+2.1 


-0.0, 


+0.0 


-0.0, 


+0.0 


Total systematic 


-6.9. 


+7.6 


-5.1. 


+4.7 


-6.2, 


+6.5 


Cross section [45] 


+ 1.1 




±1.2 




±0.5 




Total statistical 


±4.2 




±3.7 




±9.3 





TABLE XX: Systematic uncertainties (%) on fluxes for the energy- 
constrained analysis of the salt data set. Note that "const." denotes an 
energy-independent systematic component and "E dep" an energy- 
dependent part. 




FIG. 29: Flux of /j + t neutrinos versus flux of electron neutri- 
nos. CC, NC and ES flux measurements are indicated by the filled 
bands. The total 8 B solar neutrino flux predicted by the Standard So- 
lar Model [13] is shown as dashed lines, and that measured with the 
NC channel is shown as the solid band parallel to the model predic- 
tion. The narrow band parallel to the SNO ES result correponds to 
the Super- Kamiokande result in [9]. The intercepts of these bands 
with the axes represent the ±ltr uncertainties. The non-zero value 
of (pjjr provides strong evidence for neutrino flavor transformation. 
The point represents <p e from the CC flux and <^ T from the NC-CC 
difference with 68%, 95%, and 99% CL. contours included. 



Source 


NC uncert. (%) 


CC uncert. (%) 


ES uncert. (%) 


Energy scale (const.) 


-0.3, 


+0.7 


-3.7, 


+3.9 


-1.8, 


+ 1.6 


Energy scale (E dep.) 


-0.9, 


+ 1.0 


-1.0. 


+ 1.0 


-0.2, 


+0.2 


Energy radial bias 


-0.1. 


+0.1 


-2.5, 


+2.6 


-1.0, 


+0.9 


Energy resolution 


-2.1, 


+2.1 


-1.1, 


+ 1.1 


-0.6, 


+0.6 


/3\4 mean (const.) 


-2.2, 


+3.0 


-2.4, 


+2.0 


-0.5, 


+2.3 


/?i4 mean (E dep.) 


-0.2, 


+0.2 


-0.2, 


+0.2 


-0.7, 


+0.7 


/?]4 width 


-0.0, 


+0.0 


-0.1, 


+0.1 


-0.8, 


+0.8 


Radial scale (const.) 


-3.0. 


+3.3 


-2.7, 


+2.6 


-1.9, 


+2.9 


Radial scale (E dep.) 


-0.2. 


+0.2 


-1.3, 


+ 1.2 


-0.8, 


+0.8 


Vertex x 


-0.0, 


+0.1 


-0.0, 


+0.0 


-0.1, 


+0.1 


Vertex y 


-0.1, 


+0.0 


-0.0, 


+0.0 


-0.2, 


+0.2 


Vertex z 


-0.1, 


+0.1 


-0.1, 


+0.0 


-0.0, 


+0.0 


Vertex resolution 


-0.1. 


+0.1 


-0.2, 


+0.2 


-0.7, 


+0.7 


Angular resolution 


-0.2. 


+0.2 


-0.4, 


+0.4 


-4.9, 


+4.9 


Internal neutron bkgd. 


-1.9, 


+ 1.6 


-0.0, 


+0.0 


-0.0, 


+0.0 


Internal y bkgd. 


-0.2. 


+0.1 


-0.1, 


+0.0 


-0.0, 


+0.1 


Internal Cherenkov bkgd. 


-0.9, 


+0.0 


-0.8, 


+0.0 


-0.0, 


+0.0 


External Cherenkov bkgd. 


-0.2. 


+0.0 


-0.2, 


+0.0 


-0.0, 


+0.0 


Instrumental bkgd. 


-0.4, 


+0.0 


-0.3, 


+0.0 


-0.0, 


+0.0 


Neutron capture eff. 


-2.3, 


+2.1 


-0.0, 


+0.0 


-0.0, 


+0.0 


Total systematic 


-5.4, 


+5.7 


-6.2, 


+6.0 


-5.9, 


+6.6 


Cross section [45] 


±1.1 




±1.2 




±0.5 




Total Statistical 


±3.9 




±3.1 




±9.8 





Note that the uncertainties on the ratios are not normally dis- 
tributed. 

The non-Ve active neutrino component (^ T ) of the 8 B flux 
can be determined by subtracting the <p e component, as mea- 
sured by the CC flux, from the NC and ES fluxes. Whereas the 
NC measurement is equally sensitive to all active neutrinos, 
the ES measurement has reduced sensitivity to non-electron 
neutrinos in the form <t>es = 4>e + 0.1553^ T . The resulting <p^ lT 



fluxes, in units of 10 6 cm 2 s are 

^NCuncon = 3^ ± q 35 ^ +0.40 (gygt j 
^ES.uncon _ 4 35 ± j 53 ^ +0.90 (syfJ ^ 

Figure 29 shows the flux of non-electron flavor active neutri- 
nos (<p^ T ) versus the flux of electron neutrinos (<p e ). The error 
ellipses shown are the 68%, 95% and 99% joint probability 
contours for ^ T and <p e . 

Adding the constraint of an undistorted 8 B energy spectrum 
to the signal extraction yields, for comparison with earlier re- 
sults (in units of 10 5 cirr 2 s~'): 

4° c n = 1.72!^(stat)!» : ;|(syst) 
Cs" = 2.34tg(stat)t^(syst) 
Cc = 4.81^(stat)!° : 28 (syst), 

with corresponding ratios 
<A con 

-§| = 0.358 ±0.021 (stat) ™(syst) 

^NC 
.con 

^§ = 0.736 ±0.079 (stat) ™(syst), 

^ES 
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TABLE XXI: Correlation matrix for the constrained fit. External 
neutrons component is labeled EN. 





NC 


CC 


ES 


EN 


NC 


1.000 








cc 


-0.400 


1.000 






ES 


-0.073 


-0.168 


1.000 




EN 


-0.472 


-0.039 


-0.012 


1.000 



and (p^j values, in units of 10 6 cm 2 s , 

<^ T Ccon = 3.09 + 0.22 (stat) + _%™ (syst) 
^• oon = 3.97 ±1.56 (stat) (syst). 

The "energy-constrained" fit is sensitive to somewhat different 
systematic uncertainties than the unconstrained fit as shown 
in Table XX. The correlation matrix for the constrained fit, 
including correlations with the external-neutron component, 
is given in Table XXI. 

Compared to the initial salt phase results [12], some sys- 
tematic uncertainties have slightly increased. In particular, 
more detailed analyis of calibration source data during the full 
salt data set has generated larger systematic uncertainty esti- 
mates on mean flu and its energy dependence, and on angu- 
lar resolution. The combined systematic uncertainties for the 
CC and NC fluxes have not increased. The ES flux system- 
atic uncertainties have increased, but the ES measurement is 
still dominated by statistical uncertainty which has decreased 
with the increased statistics. For the (pcc/<pNC ratios, how- 
ever, the effects of flu systematic uncertainties are highly anti- 
correlated, and consequently the <PccI<Pnc systematic uncer- 
tainties are larger than that reported in [12]. 

As shown in Fig. 30 and Table XXII, these results are con- 
sistent with the pure D2O phase results [3, 4]. Comparisons 
with the 254-day salt data measured fluxes [12], also given 
in Table XXII, show some small differences for the uncon- 
strained fit case. The differences are consistent with ~ lcr sta- 
tistical fluctuations, including covariances, between the 254- 
day data set of Ref. [12] and the additional 137-day data set 
included here. Part of the difference may also be attributed to 
uncorrelated components of the systematic uncertainties be- 
tween the 254-day and 137-day data sets. 



XII. DAY-NIGHT ASYMMETRIES OF SOLAR NEUTRINO 

FLUX 

For certain ranges of mixing parameters, matter-enhanced 
neutrino oscillations predict v e regeneration inside the 
Earth [48-54]. This effect would manifest itself as a day-night 
asymmetry in the charged current and elastic scattering rates 
in SNO. For standard neutrino oscillations between active fla- 
vors, the neutral current rate should not vary between day and 
night. An observation of a day-night asymmetry in the neutral 



Phase I (306 days) 
constrained fit 

Phase II (391 days) 
constrained fit 

Phase II (391 days) 




unconstrained fit 



j 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 

1.3 1.4 1.5 1.6 1.7 1.8 1.9 

<|> cc (x 10 6 cnr 2 s 1 ) 



Phase I (306 days) 
constrained fit 



Phase II (391 days) 
constrained fit 



Phase II (391 days) 




unconstrained fit 



j 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 

3.5 4 4.5 5 5.5 6 

<Sf NC (x 10 6 cm 2 s 1 ) 



Phase I (306 days) 
constrained fit 



Phase II (391 days) 
constrained fit 



Phase II (391 days) 




unconstrained fit 





1.4 1.6 1.8 2 2.2 2.4 2.6 2.8 

<)> ES (x 10 6 cm 2 s 1 ) 

FIG. 30: Comparison of phase I (306 days) and phase II (391 days) 
flux results. For each case the inner error bar represents the statistical 
uncertainty while the full error bar represents the combined statistical 
and systematic uncertainty. 
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TABLE XXII: Constrained and unconstrained flux results from phase 
I and phase II SNO data sets in units of 10 s cirr 2 s _1 . Note that the 
phase I T cff threshold was lower than the phase II threshold. 

Constrained fit 



Data Set 



1.76 



Phase I (306 days) [4] 

Phase II (254 days)[12] 1.70*; 
Phase II (391 days) 



con /h con /A con 

CC ^NC ™ES 

+0H5+TOJ9 «: nq+U.44+0.46 9 -jq+U.24+0.12 

-0.05-0.09 



2.39 + 



+0.07+0.09 
0.07-0.10 
1 79+0.05+0.11 
i - '^-0.05-0. 11 



4,01-1+0.24+0.29 9 jo+0.29+0.15 

■ y -0.24-0.27 ' -0.28-0.08 

A oi +0.19+0.28 9 -1^+0.23+0.15 

^• 01 -0.19-0.27 Z " J -0.23-0.14 



Unconstrained Fit 



Data Set 



juncon 

''cc 



juncon 

''nc 



juncon 
^ES 



Phase I (306 days) [4] 
Phase II (254 days) [12] l.59 + _° m 

Phase II (391 days) 



6 42+1.57+0.55 

-1.57-0.58 
^ 91 +0.27+0.38 991+0.31+0.10 



0.07-0.08 -0.27-0.38 ^"^ -0.26-0.10 

1 (CO+0.06+0.08 4 Q4+0.21+0.38 9 tc+0.22+0.15 



TABLE XXIII: Event totals and rates for the day and night data sets. 



Events Rate (day ) 
"Day 2134 12.09 ± 0.26 
Night 2588 12.04 ± 0.24 



TABLE XXIV: Day-night integral fluxes and asymmetries from a 
shape-unconstrained signal extraction. Fluxes are in units of 10 6 
neutrinos crrT 2 sec -1 . The systematic uncertainties on the day and 
night fluxes include large correlated systematics that cancel in the 
day-night asymmetry ratio A. 



Signal 


Day flux 


Night flux 


A 


CC 
NC 
ES 


1.73 ±0.09 ±0.10 
4.81 ±0.31 ±0.39 
2.17 ±0.34 ±0.14 


1.64 ±0.08 ±0.09 
5.02 ±0.29 ±0.41 
2.52 ±0.32 ±0.16 


-0.056 ± 0.074 ± 0.053 
0.042 ± 0.086 ± 0.072 
0.146 ±0.198 ±0.033 



current rate could be evidence for an admixture of sterile neu- 
trinos, or for unexpected matter interactions inside the earth. 
To search for day-night asymmetries in solar neutrino reaction 
rates, the day-night asymmetry ratio A = 2(<Pn-(Pd)/(4>n+<Pd) 
is constructed using the day and night fluxes (p D and (f> N for 
each reaction. The day and night neutrino fluxes at SNO can 
be measured either allowing for a neutral current day-night 
asymmetry, or constraining A N c = as predicted by standard 
neutrino oscillations with only active flavors. 

SNO's data set was divided into day and night portions, de- 
fined by the Sun being above or below the horizon respec- 
tively. Separate day and night signal probability distributions 
were built for CC, ES, and NC events from Monte Carlo sim- 
ulations that properly included the live time exposure of the 
data set. The day-night analysis used the same event selec- 
tion, analysis cuts, and background estimates as the integral 
flux and energy spectral analyses. Each background was di- 
vided between the day and night data sets according to its 
measured diurnal rate asymmetry. Fits of the signal and back- 
ground PDFs to the data sets determined the night and day 
neutrino fluxes for each type of interaction. In addition to day 
and night live time corrections, the fluxes were corrected for 
seasonal variations in the neutrino rate due to the eccentricity 
of the Earth's orbit. The eccentricity corrections were deter- 
mined from Monte Carlo simulations. 

Because most systematics cancel when forming a day-night 
ratio, the dominant uncertainties on day-night ratios are sta- 
tistical. To avoid introducing statistical biases into the anal- 
ysis, the entire data set was divided into a 20% "open" por- 
tion and an 80% "closed" portion. The 20% open fraction 
was sampled uniformly from each run. All analysis cuts and 
procedures were developed and tested based on measurement 
of day-night ratios on only the 20% open data set. Then the 
analysis procedures were frozen, and day-night asymmetries 
were calculated for the 391 -day data set. Day-night results for 
the 20% open data set and for the full data set are statistically 
consistent. Only results for the total data set are reported here. 

Systematic uncertainties on the day-night asymmetries are 
described in Sec. VIII. The effects of day-night differences in 



detector response have been determined by perturbing the day 
and night signal PDFs by the detector response uncertainty for 
each systematic. This perturbation technique is described in 
Sec. IX. Uncertainties in backgrounds have been propagated 
by varying the amplitude and the day-night rate asymmetry of 
each background in the flux fits. 



A. Total Event Rate 

The simplest day-night analysis is to compare the total 
event rates (signals + backgrounds) between day and night. 
Table XXIII shows these results. The day and night rates are 
statistically consistent. Because the external neutron back- 
ground is determined from fits to the data itself, backgrounds 
cannot be subtracted from the raw event rates without doing a 
full signal extraction fit. 



B. Model-Independent Day-Night Asymmetries 

The most general day-night analysis is to fit for the day 
and night neutrino fluxes separately, placing no constraint on 
Anc and making no assumption about the energy dependence 
of the v e oscillation probability. The results include day and 
night NC fluxes, and separate day and night CC energy spec- 
tra. 

Table XXIV presents the day and night integral fluxes from 
the shape-unconstrained analysis. Each pair of day-night 
fluxes shares some large common systematics, as calculated 
for the integral flux analysis in Sec. XI. The day and night 
fluxes are statistically independent, however. The asymme- 
try ratio A for each flux includes a statistical uncertainty, and 
a systematic uncertainty due to day-night specific effects as 
described in Sec. VIII. All asymmetries are consistent with 
zero. 

Figure 31(a) shows the value of Acc in each energy bin 
Ace,;- Overlaid is the expectation for the previous best-fit mix- 
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FIG. 3 1 : Day-night asymmetries on each CC energy bin as a function 
of electron energy. Panel (a) shows the case in which no constraint is 
made on Anc- Panel (b) shows the case in which Anc is constrained 
to zero. Uncertainties are statistical only. The vertical lines in each 
figure show the expectation for Am 2 = 7x 1(T 5 eV 2 and tan 2 6 = 0.40. 



ing parameters Am 2 = 7 x 1CT 5 eV 2 and tan 2 6» = 0.40 [12]. 
The dependence of Aqc on CC electron energy is consistent 
with this expectation, but is also consistent with no day-night 
asymmetries. 

Although the day and night fluxes are statistically indepen- 
dent, the CC, ES, and NC fluxes for either day or night are sta- 
tistically correlated since they are produced from a common 
fit. As a result, Acc, ^nc, and Aes are statistically correlated, 
with correlation coefficients given by: 

p(CC, NC) = -0.532 
p(CC,ES) = -0.147 
p(ES, NC) = -0.064. 

Similarly, the Aqcj are modestly correlated between differ- 
ent energy bins at lower energies due to their common covari- 
ance with Anc- Figure 32 is a contour plot of Anc versus Acc, 
illustrating the covariances. 

Table XXV lists each systematic uncertainty on the inte- 
gral flux asymmetries from the shape-unconstrained day-night 
analysis. The largest systematics are uncertainties on diurnal 



TABLE XXV: Systematic uncertainties on day-night asymmetries 
for the shape-unconstrained signal extraction. For presentation, un- 
certainties have been symmetrized and rounded. 



Systematic 


Aqc unceit. 


A^c unceit. 


Aes uncert. 


Diurnal energy scale 


0.004 


0.015 


0.007 


Directional energy scale 


0.001 


0.000 


0.014 


Long-term energy scale variation 


0.002 


0.010 


0.001 


Diurnal energy resolution 


0.003 


0.006 


0.004 


Directional energy resolution 


0.001 


0.001 


0.003 


Diurnal vertex shift 


0.008 


0.012 


0.007 


Directional vertex shift 


0.000 


0.000 


0.003 


Diurnal vertex resolution 


0.002 


0.006 


0.002 


Directional vertex resolution 


0.000 


0.000 


0.001 


Diurnal isotropy 


0.050 


0.064 


0.017 


Directional isotropy 


0.002 


0.002 


0.004 


Long-term isotropy variation 


0.014 


0.015 


0.006 


Directional angular resolution 


0.001 


0.001 


0.020 


Live time 


0.000 


0.000 


0.000 


Cut acceptance 


0.003 


0.004 


0.003 


External Cherenkov tail bkgd. amplitude 


0.002 


0.003 


0.000 


External Cherenkov tail bkgd. asymmetry 


0.003 


0.004 


0.000 


Internal Cherenkov tail bkgd. amplitude 


0.001 


0.000 


0.000 


Internal Cherenkov tail bkgd. asymmetry 


0.001 


0.001 


0.000 


Internal neutron bkgd. amplitude 


0.000 


0.003 


0.000 


Internal neutron bkgd. asymmetry 


0.000 


0.015 


0.000 


Internal y bkgd. amplitude 


0.000 


0.000 


0.000 


Internal y bkgd. asymmetry 


0.000 


0.000 


0.000 


Isotropic AV bkgd. amplitude 


0.001 


0.001 


0.000 


Isotropic AV bkgd. asymmetry 


0.002 


0.002 


0.000 


Instrumental bkgd. amplitude 


0.000 


0.001 


0.000 


Instrumental bkgd. asymmetry 


0.001 


0.002 


0.000 


Total 


0.053 


0.072 


0.033 
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TABLE XXVI: Day-night integral fluxes from a shape-unconstrained 
signal extraction, with the constraint Anc = 0. Fluxes are in units of 
10 6 neutrinos citT 2 sec -1 . 



> 0.25 _ 



Signal 


Day flux Night flux 


Asymmetry 


CC 


1.71 ± 0.08 ± 0.09 1 .65 ± 0.08 ± 0.09 


-0.037 ± 0.063 ± 0.032 


ES 


2.18 ±0.34 ±0.14 2.53 ±0.32 ±0.16 


0.153 ±0.198 ±0.030 


NC 


4.93 ±0.21 ±0.36 


A NC = 



variations of the isotropy parameter f3\4, and diurnal variations 
in energy scale and vertex shift. For the ES rate, directional 
systematics are significant. However, the overall uncertainties 
on all asymmetries are ultimately limited by statistics. 



C. Shape-Constrained Day-Night Asymmetries 

A variant of the preceding analysis is to constrain the day 
and night v e energy spectra to follow an undistorted 8 B shape. 
This corresponds to an energy-independent oscillation prob- 
ability that varies between night and day. The NC rate was 
again allowed to vary in the fit. It should be noted that stan- 
dard neutrino oscillations with mixing parameters in the LMA 
region do not predict energy-independent day-night asymme- 
tries. The derived asymmetries under the assumption that the 
CC and ES energy spectra are undistorted, but allowing them 
to have different normalizations between night and day, are 

A cc = -0.021 +0.063 ±0.035 

Anc = 0.018 ±0.079 ±0.052 

A ES = 0.066 ±0.198 ±0.057. (17) 



D. Day-Night Asymmetries with the Constraint A NC = 



o o 
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FIG. 33: (a) Day and night extracted CC energy spectra with statisti- 
cal uncertainties only. Bin values are expressed in units of equivalent 
8 B flux, normalized such that the sum of the flux bin values above 5.5 
MeV equals the total integral 8 B neutrino flux above MeV, as deter- 
mined for the day and night integral fluxes quoted in section XII(D) 
(see Appendix A), (b) Difference night - day between the spectra. In 
both figures, the final bin extends to 20 MeV. 



In the previous two subsections the NC flux was allowed 
to vary in the fit between the day and night data sets. Under 
the standard picture of matter-enhanced neutrino oscillations, 
Anc should be zero. This prediction has been confirmed by 
the results of the previous two subsections. When determin- 
ing the best estimate of the day-night asymmetry on the elec- 
tron neutrino flux, assuming standard neutrino oscillations, it 
is appropriate to constrain Anc- This constraint has been ap- 
plied by simultaneously fitting the day and night data sets, not 
allowing 0nc to vary between night and day. 

The additional constraint of Anc — reduces the statisti- 
cal uncertainties on Acc and Aes- It also produces a modest 
covariance between the day and night fluxes, due to their com- 
mon covariance with 0nc- In contrast, without a constraint on 
Anc the day and night neutrino fits are statistically indepen- 
dent. 

The day and night neutrino fluxes were fit in a shape- 
unconstrained analysis, requiring Anc = 0. Table XXVI gives 
the day and night integral fluxes from this fit, and the NC flux. 
No statistically significant asymmetries are observed. Forcing 
Anc — results in some reduction in |Accl, as expected from 



the anti-correlation of CC and NC event totals in the signal 
extraction. 

Figures 31(b) and 33 show the CC asymmetry as a function 
of electron energy and the day and night CC energy spectra, 
binned by electron energies. The additional constraint does 
not significantly change the results. Table XXVII presents the 
systematic uncertainties on Acc and A E s for this analysis. 



E. Shape-Constrained Day-Night Asymmetries with the 
Constraint A NC — 

For the sake of completeness the analysis of Sec. XII C has 
been repeated with the additional constraint that Anc = 0. The 
results are 

A cc = -0.015 ±0.058 ±0.027 

A ES = 0.070 ±0.197 ±0.054. (18) 
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TABLE XXVII: Systematic uncertainties on day-night asymmetries 
for the shape-unconstrained signal extraction, with the constraint 
A NC = 0. For presentation, uncertainties have been symmetrized 
and rounded. 



Systematic 


Acc uncert. 


Aes uncert. 


Diurnal energy scale 


0.009 


0.009 


Directional energy scale 


0.001 


0.014 


Long-term energy scale variation 


0.006 


0.002 


Diurnal energy resolution 


0.002 


0.004 


Directional energy resolution 


0.001 


0.002 


Diurnal vertex shift 


0.013 


0.009 


Directional vertex shift 


0.000 


0.003 


Diurnal vertex resolution 


0.001 


0.002 


Directional vertex resolution 


0.000 


0.001 


Diurnal isotropy 


0.022 


0.009 


Directional isotropy 


0.001 


0.005 


Long-term isotropy variation 


0.013 


0.003 


Directional angular resolution 


0.002 


0.019 


Live time 


0.000 


0.000 


Cut acceptance 


0.003 


0.003 


External Cherenkov tail bkgd. amplitude 


0.002 


0.000 


External Cherenkov tail bkgd. asymmetry 


0.003 


0.000 


Internal Cherenkov tail bkgd. amplitude 


0.001 


0.000 


Internal Cherenkov tail bkgd. asymmetry 


0.001 


0.000 


Internal neutron bkgd. amplitude 


0.002 


0.001 


Internal neutron bkgd. asymmetry 


0.007 


0.002 


Internal y bkgd. amplitude 


0.000 


0.000 


Internal y bkgd. asymmetry 


0.000 


0.000 


Isotropic AV bkgd. amplitude 


0.001 


0.000 


Isotropic AV bkgd. asymmetry 


0.002 


0.000 


Instrumental bkgd. amplitude 


0.001 


0.000 


Instrumental bkgd. asymmetry 


0.001 


0.000 


Total 


0.032 


0.030 



F. Combined Day-Night Asymmetries from the Salt and D 2 
Data Sets 

During the first phase of the SNO experiment, the asym- 
metry on the electron neutrino flux was measured to be A e = 
0.070 + 0.049+[}^, assuming a standard 8 B shape and con- 
straining Anc = 0. Although an asymmetry ratio formed from 
two normally-distributed variables is not necessarily normally 
distributed, for the case of SNO's day and night fluxes a nor- 
mal distribution is an excellent approximation for the true dis- 
tribution we calculate for Acc or A e . The asymmetry results 
from SNO's first phase are statistically independent of the re- 
sults from the salt data set and statistical uncertainties dom- 
inate over systematics for the asymmetries. Hence, combin- 
ing values and uncertainties for Acc from Sec. XII E and A e 
from above can be done trivially to produce a combined, al- 
beit model-dependent, day-night asymmetry of A eiCom b me d = 
0.037 + 0.040. A future SNO analysis will address the issue 
of doing a joint shape-unconstrained fit to the D2O and salt 
data sets. 

The Super-Kamiokande collaboration has measured the 
day-night rate asymmetry of ES interactions above 5 MeV 
to be A ES = 0.021 + 0.020+°;^ [55]. Because ES inter- 
actions can be initiated by either v e or v^, the day-night 
asymmetry for ES events is diluted by a factor of (<f> e + 
O.1576(0 fo , - (p e ))/((l - 0.1576)<7> e ) = 1.55. Assuming an 
energy-independent conversion mechanism and only active 
neutrinos, the Super-Kamiokande result scales to a <p e asym- 
metry of A e sK = 0.033 ± 0.031^;^. Combining the SNO 



D2O and SNO salt values for A e with the equivalent Super- 
Kamiokande value A Ci sk gives A eiCom bi ne d = 0.035 ± 0.027. 
This result explicitly assumes a shape-constrained 8 B spec- 
trum for SNO and ignores energy dependence of the oscilla- 
tion probability over the energy range in question. 



XIII. MSW INTERPRETATION OF RESULTS 

The observation of a substantially suppressed v e flux with 
the CC reaction in SNO compared to the total active flux mea- 
sured by the NC reaction in SNO provides clear evidence for 
neutrino flavor change that can be analyzed in terms of neu- 
trino oscillations. Constraints on neutrino mixing parameters 
can be derived by comparing neutrino oscillation model pre- 
dictions with experimental data, as has been done in, for ex- 
ample, [53, 56-60] and in previous SNO analyses [5, 12]. 

A two-flavor, active neutrino oscillation model has two 
parameters: Am 2 , the difference between the square of the 
masses of the relevant eigenstates of propagation for the neu- 
trinos, and tan 2 9 that quantifies the strength of the mixing be- 
tween flavor and mass eigenstates. Note that the three-flavor 
mixing matrix element U e 2 can be written as cos #13 sin #12 
[15], which is approximately equal to sin 6 for two-flavor so- 
lar neutrino oscillations when #13 is small and when Am 2 ol <s 
Am 2 m . The MSW effect [16] can result in neutrinos above a 
few MeV emerging from the Sun essentially as a pure V2 state 
(e.g. for oscillation parameters in the Large Mixing Angle re- 
gion). To the degree to which this statement is true, SNO's 
</>cc/<Pnc ratio, a direct measure of the v e survival probability, 
is also a direct measure of \U e 2\ 2 and thus should be approxi- 
mately equal to sin 2 6. For the sake of comparison with other 
past and present oscillation analyses, this present work still 
employs tan 2 8 to quantify the mixing angle for solar neutrino 
oscillations. 

For each pair of parameters, the oscillation model predicts 
the expected rates in the CI [6] and Ga experiments [7, 8], 
Super-Kamiokande zenith spectra [9], and SNO rates and 
spectra. The model prediction accounts for MSW propaga- 
tion of neutrino states through dense matter in the Sun and the 
Earth and so allows for the regeneration of v e flavor for neutri- 
nos passing through the Earth at night. A global^ 2 calculation 
can be performed; best-fit parameters can be determined and 
allowed parameter regions can be identified using A^ 2 confi- 
dence levels for two degrees of freedom. The same neutrino 
oscillation model can also predict rates and spectra for the 
KamLAND experiment [17, 61], assuming CPT invariance 
(since KamLAND detects v e ). The likelihood values from a 
KamLAND oscillation analysis can be easily combined with 
that from the global solar neutrino analysis to further restrict 
the allowed oscillation parameter space. 

For the analysis presented in this paper, earlier data from 
SNO-I (pure D 2 phase) have been included. SNO-I day and 
night spectra have been interpreted in a similar manner as be- 
fore [5, 12, 62]. Summed spectra (CC+ES+NC+background) 
predictions were compared to the number of counts in each 
spectral bin from the SNO-I data, for both day and night. A 
minor improvement in the analysis of this earlier data is that 
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this part of the calculation now includes energy-dependent v-d 
radiative corrections for the CC reaction. Previously the CC 
radiative correction was included as an energy-independent 
factor. 

In [12], salt phase "fluxes" (i.e., CC, NC and ES fluxes in- 
ferred from rates) were added to the global^ 2 analysis. The 
present work has a new oscillation analysis using data from 
the 391 -day data set of the salt phase which have been ana- 
lyzed and extracted as CC spectra and NC and ES integrated 
fluxes, separately for day and night. This information is in- 
cluded in the global^ 2 analysis in lieu of just salt phase fluxes 
inferred from rates. This allows CC spectral shape informa- 
tion and day-night rate asymmetry information from the salt 
phase to be included in the global oscillation analysis. CC- 
NC separation is preserved in this analysis since the SNO un- 
constrained signal extraction utilized information from event 
isotropy fiw and angular correlation cos 9 Q distributions for 
separating the salt NC and ES fluxes from the CC spectra. 

SNO's unconstrained signal extraction produced two 
19x19 statistical covariance matrices (one for day and one for 
night) for 17 spectral bins of the CC spectrum, starting from 
5.5 MeV kinetic energy up to 13.5 MeV, in 0.5 MeV steps, 
with one extra bin integrating from 13.5-20.0 MeV, plus the 
NC and ES fluxes. These statistical covariance matrices are 
required in the calculation ofx 2 and are available in Appendix 
A. Day and night data are statistically independent from each 
other and the results with no constraint on Anc were used. 

Systematic uncertainties also have bin-to-bin correlations 
and unlike the statistical correlations from SNO's signal ex- 
traction may also include correlations that extend across day 
and night spectra. Experimental spectral shape systematic un- 
certainties were described in Section X and were included in 
this oscillation analysis. The uncertainty in the shape of the 
8 B neutrino spectrum has also been included in this x 1 analy- 
sis. The 8 B spectrum used in our model is the one from [38]; 
however, the more generous uncertainties from [63] were em- 
ployed in the systematics calculation in our^- 2 analysis. 

Day-night systematics, though small, were also included 
in the global x 1 analysis. The significant day-night system- 
atics are diurnal energy scale, long-term energy scale varia- 
tion, diurnal vertex shift, diurnal isotropy variation, long-term 
isotropy variation, and internal neutron background asymme- 
try. Other day-night systematics discussed in this paper are 
smaller in magnitude and were averaged together in the^ 2 cal- 
culation. Note that some directional systematics have a non- 
negligible effect on the day-night asymmetry of ES events; 
however, the impact of the day-night asymmetry of ES events 
in SNO on the oscillation analysis is not that significant so 
combining these systematics is also reasonable. The tech- 
nique for including systematic uncertainties and bin-to-bin 
correlations in the x 2 analysis is the conventional one, as in 
[64]. Thus^ NO u from SNO's 391 -day data set is defined as: 
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2 \ 1 /T/data T/modeKr 2 -1 /T/data T/modeK /in\ 

*SNO-II = 2j (F < ~ Y i ^i/tOt)] ( Y j ~ Y j )' ( 19 > 

U=l 

where T^ ata is the SNO experimental value in one of the 17 CC 
spectral bins, or the NC or ES flux, day or night, and y™ odel is 



the model predicted value for bin i based on the neutrino os- 
cillation hypothesis and the set of parameters being evaluated. 

The error matrix for the calculation <x 2 .(tot), is composed of 
statistical and systematic components: 

cr 2 (tot) = (T 2 (stat) + (T 2 (syst), (20) 

with cr 2 .(stat) containing the elements from the statistical co- 
variance matrices from SNO's unconstrained signal extraction 
and <x 2 .(syst) containing contributions from systematic uncer- 
tainties. The spectral systematics error matrix is formed from 
the partial derivatives that relate the rate y model in the rth bin to 
the uncertainty in each one of the K spectral systematics S%: 

, dYi dYj , 

cr 2 (syst)=y — i-/(A5,) 2 , (21) 

where AS * is the uncertainty estimated for spectral systematic 
St. Note that all systematic uncertainties have an effect on 
the extracted CC spectra, and possibly an energy-dependent 
effect; thus, all systematics are spectral systematics. In this 
standard^ 2 treatment, bin-to-bin correlations are included for 
the systematics; however, possible correlations among the var- 
ious systematic uncertainties were neglected. 




tan 2 9 

FIG. 34: SNO-only neutrino oscillation analysis, including pure D2O 
phase day and night spectra, and salt extracted CC spectra, NC and 
ES fluxes, day and night. The 8 B flux was free in the fit; hep solar 
neutrinos were fixed at 9.3 x 10 3 cnT 2 s~' . The star is plotted at the 
best-fit parameters from the^ 2 analysis, listed in Table XXVIII. 

Figure 34 shows the allowed regions for neutrino oscillation 
parameters when only SNO data (SNO-I and SNO-II) are an- 
alyzed. The inclusion of CC spectral data, improved measure- 
ment of the NC flux from the larger data set, and the addition 
of separate day and night results compared with [12] produce 
slightly smaller allowed ranges of parameters. The best-fit 
parameters from a SNO-only analysis are: Am 2 = 5.0 x 10~ 5 
eV 2 , tan 2 6 = 0.45, fs — 5.1 1 x 10 6 cm" 2 s" 1 , which is the total 
active 8 B solar neutrino flux, a free parameter during^ 2 min- 
imization. The best-fit^ 2 is 68.9 for 69 degrees of freedom in 
the SNO-only oscillation analysis. 
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FIG. 35: (a) Global neutrino oscillation analysis using only so- 
lar neutrino data, and (b) including KamLAND 766 ton-year data. 
The solar neutrino data included SNO's pure D 2 phase day and 
night spectra, SNO's salt phase extracted day and night CC spectra 
and ES and NC fluxes, the rate measurements from the CI, SAGE, 
Gallex/GNO, and SK-I zenith spectra. The 8 B flux was free in the 
fit; hep solar neutrinos were fixed at 9.3 x 10 3 cirT 2 s _I , The stars 
are plotted at the best-fit parameters from the x 2 analysis, listed in 
Table XXVIII. 



but this shift is perfectly consistent with the global solar neu- 
trino constraints and gives ax 2 = 1 13.6 for the solar neutrino 
part of the calculation. A summary of the best-fit oscillation 
parameters and their ranges within the allowed LMA regions 
appears in Table XXVIII. SNO data are providing strong con- 
straints on the mixing angle. 



TABLE XXVIII: Best-fit neutrino oscillation parameters. Uncertain- 
ties listed are ±lcr for the 2-D parameter regions (and only within 
the LMA region for the SNO-only analysis). 



Oscillation analysis 


Am 1 (lfr 5 eV 2 ) 


tan 2 


SNO-only 

Global solar 

Solar plus KamLAND 


5 +6 - 2 
6.5^1 
8.0$ 


45+0.11 
u -^ J -o.io 
45+ 09 

" -0.08 

45+ 09 

"'^ -0.07 



Compared to [12] the inclusion of the 391 -day salt data 
set (with spectral and day-night information) in the oscilla- 
tion analysis moves the allowed oscillation region to slightly 
larger mixing angles. This is due to the larger central value 
of the <pcc/4>NC ra tio found in the present analysis. The 2004 
KamLAND data [61] have already tightly constrained the pa- 
rameter Am 2 . In terms of individual uncertainties the results 
become Am 2 = 8.0+^ x 10~ 5 eV 2 and G = 33.9^ degrees, 
where the uncertainties were obtained as 1 -dimensional pro- 
jections of the respective parameter while marginalizing the 
uncertainties in the other. 



TABLE XXIX: Comparison of SNO total active 8 B solar neutrino 
flux measurements and solar model predictions. 



Source 


Total S B Flux (10 6 cr' s" 1 ) 


SNO pure D 2 phase NC 


5.09«"(stat.)«;«(syst.) 


above, energy unconstrained 6.42 ± 1.57(stat.)+?~(syst.) 


SNO salt phase NC 


4.94 ± 0.2 l(stat.) + °-^(syst.) 


SNO salt day NC 


4.81 ± 0.31(stat.) ± 0.39(syst.) 


SNO salt night NC 


5.02 ± 0.29(stat.) ± 0.41(syst.) 


SNO-only oscillation lit 


5.11 


global solar fit 


5.06 


solar plus KamLAND fit 


4.93 


BS05(OP) [13] 


5.69 ±0.91 


BS05(AGS,OP) [13] 


4.51 ±0.72 


BP04 [11] 


5.79 ± 1.33 


BP2000 [10] 


5 05 +lm 


TC04 tac A [14] 


4.25 


TC04 seismic [14] 


5.31 ±0.6 



The top panel in Fig. 35 shows the allowed region for 
a global oscillation analysis that included data from all so- 
lar neutrino experiments. The best-fit oscillation parameters, 
with lcr uncertainties on the 2-dimensional parameter region 
given, are Am 2 = 6.5+^ x 10~ 5 eV 2 , tan 2 (9 = 0.45^, with 
a best-fit^ 2 = 113.1 for 116 degrees of freedom in the global 
solar x 1 analysis. The lower panel shows the results of the 
analysis when the 766 ton-year data from KamLAND [61] 
were also included. The best-fit parameters from the global 
solar plus KamLAND analysis are: Am 2 = 8.O+JJ 4X 10~ 5 eV 2 , 
6 = 33.9! 2 !} degrees, f B = 4.93 x 10 6 cirT 2 s -1 , where the la- 
uncertainties on the 2-dimensional parameter region are given. 
The inclusion of KamLAND data shifts the best-fit Am 2 value 



The total active B solar neutrino flux, measured by the NC 
reaction, has been presented in several ways in SNO analyses. 
Table XXIX lists SNO measured (or fit) values and fluxes pre- 
dicted by solar models. In the first row, the SNO NC flux was 
extracted assuming an undistorted 8 B spectrum (for the null 
hypothesis test). All subsequent values in the table are free 
from that assumption. The salt phase NC value (this work) 
is the most precise and appropriate one to compare with solar 
models. The agreement between solar models and this mea- 
surement is good. 

Based on the best-fit parameters from the global solar plus 
KamLAND analysis, the predicted CC electron energy spec- 
trum is determined. In Fig. 36, this prediction is compared 
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to the measured CC spectrum. The x 1 between the extracted 
spectrum and the expected shape for the best-fit LMA param- 
eters, calculated with all statistical correlations and system- 
atic uncertainties as described above for the global oscillation 
analysis, is 27.2 for 16 degrees of freedom (17 spectral bins 
minus a floating normalization factor). The probability of ob- 
serving a x 2 > 27.2 under the assumption that the data are 
drawn from the expected LMA spectrum is 3.9%. 



Data 

Systematic uncertainties 
SSM 8 B model shape 
LMA 8 B model shape 




sented. The salt additive enables a statistical separation of NC 
events from CC and ES events by measuring event isotropy. 
In addition to new results for integral fluxes, energy spec- 
tral information from the CC reaction is presented, with com- 
plete statistical and systematic uncertainties. Separate day and 
night spectra and day-night integrated-flux asymmetries are 
also presented. The flux measurements are in agreement with, 
and slightly more precise than, previous measurements [12]. 
The energy spectrum derived from the CC reaction is consis- 
tent with the expected spectrum assuming an undistorted 8 B 
shape and also with the predicted spectrum corresponding to 
the best-fit LMA parameters for a global oscillation analysis 
with solar neutrino and KamLAND reactor neutrino data in- 
cluded. Within uncertainties, no significant day-night asym- 
metries are observed as expected for the best-fit LMA solu- 
tion. Detailed MSW fits find a single allowed region in oscil- 
lation parameter space and tightly constrained values for Am 2 
and the mixing angle 9. These data provide further confirma- 
tion of flavor change for solar neutrinos and for the oscillation 
of massive neutrinos as the dominant flavor change mecha- 
nism. The total flux of all active neutrino types for 8 B solar 
neutrinos is in agreement with the most recent solar model 
calculations [13, 14]. 



FIG. 36: Extracted CC T cS spectrum compared to that predicted with 
the best-fit LMA parameters. Only statistical uncertainties are shown 
in the data spectrum. The band on the undistorted 8 B model shape 
represents the lcr uncertainty determined from detector systematic 
uncertainties. The predicted spectrum is normalized to the same 
number of counts as the data spectrum. Note that the data points, 
especially the first three points, are statistically correlated as well as 
having correlated systematics as indicated by the error band. 



XIV. SUMMARY 

An extensive analysis of the data from the full running pe- 
riod with salt added to the heavy water in SNO has been pre- 
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APPENDIX A: SNO DATA IN AN OSCILLATION ANALYSIS 

The following tables contain a subset of SNO's salt phase 
results and some supporting information that may be needed 
in a neutrino oscillation analysis. The values in these tables 
include outputs from the signal extraction (described earlier 
in this paper) and results from systematic uncertainty stud- 
ies. Day and night values were extracted separately without 
any constraints on the NC rate asymmetry or energy spectrum 
shapes. 

Table XXX contains day and night CC spectra, expressed 
as fluxes. The meaning of these fluxes (for example the day 
flux in the 6.0-6.5 MeV bin of 0.182 x 10 6 cnr 2 s _1 ) is that 
the number of events SNO observed in the salt day data set 
attributed to CC interactions by the signal extraction, with an 



electron kinetic energy between 6.0-6.5 MeV, is equal to the 
number of all CC events that would be observed above kinetic 
energy 5.5 MeV, if the integral flux (from zero to endpoint) of 
v e had the value of 0.182xl0 6 crrT 2 s _1 and had an undistorted 
8 B spectral shape [38]. The 8 B spectral shape aspect of this 
definition is only for normalization. There is no assumption 
of any spectral shape when extracting the actual number of 
events in each bin for SNO's salt phase. This normalization 
was chosen so that the sum of the values for all bins equals the 
integral day and night 8 B fluxes quoted in Section XII B. 

When calculating the theoretical CC flux and spectra for a 
set of oscillation parameters, for comparison with SNO data, 
one should be aware of the above definition. Thus, a model 
prediction for the value in the above mentioned example bin 
would be: 



J 



Jo 



Jo° I°° £ ${E v )Pee{E v )§-{E v , T e )R(T e , r eff ) dE v dT e dT eff 
' 1°° T E <KE v )§;(Ey, T e )R(T e , T eB ) dE v dT e dT eff 

I 



(f>(E v ) dE v 



(Al) 



where P ee is the survival probability for a v e produced in the 
Sun to be detected as a v e , <f>(E v ) is the flux of 8 B solar neu- 
trinos as a function of neutrino energy, jf- is the differential 
cross section for the CC reaction and R(T e , r eff ) is the energy 
response function, 



R(T e , T eS ) = 



1 



exp 



TMTj 



2a\ 



(A2) 



and T e is the true recoil electron kinetic energy and T e ff is the 
observed electron kinetic energy, with resolution, 



cr T (T e ) = -0.131 + 0.383 y/F e + 0.03731 T e 



(A3) 



in units of MeV. 

Tables XXXII and XXXIII contain the statistical correla- 
tion coefficients from SNO signal extraction, for salt phase 
day and night data respectively. The numbering of the CC 
spectral bins listed in these tables follows the ordering of en- 
ergy bins as listed in Table XXX. These correlation coeffi- 
cients are necessary to include in an oscillation analysis that 
includes SNO CC spectral information. 



The systematic uncertainties for the SNO extracted CC 
spectrum are listed in Table XXXIV. These uncertainties (in 
percent) can be used as the partial derivatives in a bin-to-bin 
correlated systematics part of a x 2 calculation, as described 
in Section XIII. The spectral shape of these uncertainties can 
also be used as the shape of the related day-night asymmetry 
systematic. The magnitude of a day-night spectral systematic 
can be estimated by taking the ratio of the size of the day-night 
flux asymmetry uncertainty and the total flux uncertainty, for 
each relevant systematic, and using this ratio in each bin to 
scale that spectral systematic. Note that NC-only systematics 
(internal neutron background and neutron capture efficiency 
uncertainties) are not listed in Table XXXIV since they do not 
affect the CC spectrum, but are necessary to include in an os- 
cillation analysis and can be found in Table XIX. 

Table XXXI contains the live time distribution for SNO's 
391 -day salt data set as a function of cosine of the zenith an- 
gle of the Sun. For detailed calculations of neutrino survival 
probabilities including propagation through the Earth, the live 
times at different zenith angles can be used as weighting fac- 
tors. 
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TABLE XXX: Charged-current reaction recoil electron kinetic en- 
ergy spectra from the 391 -day SNO salt phase, expressed in units of 
equivalent 8 B fluxes. The normalization is such that the sum over all 
bins equals the day or night integral 8 B solar neutrino flux above 
MeV, as determined and quoted in Section XII B. Day and night 
extracted values are listed with their statistical uncertainties from 
SNO's signal extraction. 



CC electron kinetic salt phase day salt phase night 
energy bin (MeV) ( 1 6 cm" 2 s~ 1 ) ( 1 6 cm" 2 s~ 1 ) 



5.5-6.0 


0.205 ± 0.032 


0.145 


± 0.027 


6.0-6.5 


0.182 ±0.030 


0.164 


± 0.027 


6.5-7.0 


0.153 ±0.028 


0.190 


± 0.026 


7.0-7.5 


0.226 ± 0.028 


0.180 


± 0.024 


7.5-8.0 


0.198 ±0.025 


0.178 


± 0.022 


8.0-8.5 


0.184 ±0.023 


0.164 


±0.019 


8.5-9.0 


0.124 ±0.018 


0.114 


±0.015 


9.0-9.5 


0.099 ±0.015 


0.126 


±0.015 


9.5-10.0 


0.110 ±0.015 


0.124 


±0.014 


10.0-10.5 


0.058 ±0.011 


0.067 


±0.010 


10.5-11.0 


0.070 ±0.012 


0.073 


±0.011 


11.0-11.5 


0.048 ±0.010 


0.039 


± 0.007 


11.5-12.0 


0.042 ± 0.008 


0.029 


± 0.007 


12.0-12.5 


0.0088 ± 0.0038 


0.018 


± 0.005 


12.5-13.0 


0.0082 ± 0.0040 0.015 


± 0.005 


13.0-13.5 


0.0025 ± 0.0028 


0.0042 ± 0.0025 


13.5-20.0 


0.014 ± 0.005 


0.012 


± 0.004 



TABLE XXXI: Zenith-live time distribution for the 391 -day SNO 
salt phase data set. There are sixty equally-spaced bins in cos 6 Z with 
0, being the zenith angle of the Sun (e.g., cos ft = -1 would be the 
value if the Sun was directly below the SNO detector). The amount 
of live time in each bin is listed in seconds. The bins start from 
cos ft = -1 to -0.9667 at the top left of the table and increase going 
down the column. Continuing at the top of the second column, the 
top bin has cos 9 Z = -0.5 to -0.4667. The first two columns contain 
night live times while the final two columns are day live times. 



cos ft 

-1.0 to -0.5 -0.5 to 0.0 0.0 to 0.5 0.5 to 1.0 






705975 


606806 


493795 





743962 


608735 


471616 


432202 


761182 


606873 


447927 


545975 


814881 


607682 


445073 


592399 


1107970 


616573 


473089 


566033 


887601 


617365 


491700 


570422 


777241 


625955 


453238 


539112 


728515 


651235 


458976 


527208 


698101 


682766 


466092 


571086 


672856 


752869 


449265 


594019 


652494 


844742 


433395 


613845 


639256 


633145 


388327 


622530 


624127 


576052 


255268 


672151 


615141 


570776 





684213 


608656 


521223 






TABLE XXXII: Statistical correlation coefficients from SNO's signal extraction for the salt phase day data. The numbering of the CC spectral bins follows the ordering of energy bins 
Table XXX. 







NC 




CC1 




CC2 




CC3 




CC4 




CC5 




CC6 




CC7 




CC8 




CC9 


CC10 


ecu 


CC12 


CC13 


CC14 


CC15 


CC16 


CC17 




ES 


NC 


1. 


.0000 


-0. 


.3478 


-0 


.3320 


-0 


.3281 


-0. 


.2488 


-0. 


.2000 


-0 


.1555 


-0 


.0877 


-0 


.0618 


-0. 


.0329 


-0.0152 


-0.0063 


-0.0025 


0.0003 


-0.0000 


-0.0000 


-0.0000 


-0.0000 


-0. 


0610 


CC1 


-0 


.3478 


1 


.0000 





.1528 





.1440 


0. 


1124 


0. 


.0865 





.0676 





.0425 





.0272 





.0135 


0.0061 


0.0030 


0.0011 


0.0001 


0.0000 


-0.0000 


0.0000 


0.0000 


-0, 


.0695 


CC2 


-0. 


.3320 





.1528 


1 


.0000 





.1459 


0. 


.1147 


0. 


.0873 





.0684 





.0440 





.0275 


0, 


,0135 


0.0061 


0.0031 


0.0011 


0.0001 


0.0000 


-0.0000 


0.0000 


0.0000 


-0 


.0507 


CC3 


-0 


.3281 





.1440 





.1459 


1 


.0000 


0. 


1073 


0. 


.0824 





.0644 





.0407 





.0259 





.0129 


0.0058 


0.0029 


0.0010 


0.0001 


0.0000 


-0.0000 


0.0000 


0.0000 


-0. 


.0525 


CC4 


-0. 


.2488 





.1124 





.1147 





.1073 


1. 


.0000 


0. 


.0643 





.0503 





,0322 





.0202 


0, 


,0100 


0.0045 


0.0023 


0.0008 


0.0001 


-0.0000 


-0.0000 


0.0000 


-0.0000 


-0, 


.0556 


CC5 


-0. 


.2000 


0. 


.0865 





.0873 





.0824 


0. 


.0643 


1. 


0000 





.0387 





,0243 





.0155 


0, 


,0078 


0.0035 


0.0017 


0.0006 


0.0000 


-0.0000 


0.0000 


0.0000 


0.0000 


-0. 


.0637 


CC6 


-0. 


.1555 


0. 


0676 





.0684 





.0644 


0. 


.0503 


0. 


.0387 


1 


.0000 





.0190 





.0122 





0061 


0.0027 


0.0014 


0.0005 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


-0, 


.0648 


CC7 


-0 


.0877 


0. 


.0425 





.0440 





.0407 


0. 


.0322 


0. 


.0243 





.0190 


1 


.0000 





.0077 





,0037 


0.0017 


0.0009 


0.0003 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


-0 


0595 


CC8 


-0 


0618 


0. 


.0272 





.0275 





.0259 


0. 


.0202 


0. 


0155 





.0122 





.0077 


1 


.0000 





,0024 


0.0011 


0.0005 


0.0002 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


-0 


.0430 


CC9 


-0 


.0329 





0135 





.0135 





.0129 


0. 


0100 


0. 


.0078 





.0061 





.0037 





,0024 


1 


,0000 


0.0006 


0.0003 


0.0001 


0.0000 


0.0000 


0.0000 


0.0000 


-0.0000 


-0 


.0495 


CC10 


-0 


0152 





0061 





.0061 





.0058 


0. 


.0045 


0. 


.0035 





.0027 





.0017 





.0011 


0, 


,0006 


1.0000 


0.0001 


0.0000 


0.0000 


0.0000 


-0.0000 


-0.0000 


0.0000 


-0, 


.0466 


ecu 


-0. 


.0063 


0. 


.0030 





.0031 





.0029 





.0023 


0. 


0017 





.0014 





.0009 





.0005 


0. 


.0003 


0.0001 


1.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


-0 


.0487 


CC12 


-0. 


.0025 


0. 


.0011 





.0011 





.0010 


0. 


0008 


0. 


0006 





.0005 





.0003 





.0002 


0. 


.0001 


0.0000 


0.0000 


1.0000 


0.0000 


0.0000 


0.0000 


-0.0000 


0.0000 


-0, 


.0570 


CC13 


0. 


.0003 


0. 


.0001 





.0001 





.0001 


0. 


0001 


0. 


0000 





.0000 





.0000 





,0000 


0. 


,0000 


0.0000 


0.0000 


0.0000 


1.0000 


0.0000 


-0.0000 


-0.0000 


-0.0000 


0, 


.0000 


CC14 


-0. 


.0000 


0. 


.0000 





.0000 





.0000 


-0. 


.0000 


-0. 


0000 





.0000 





,0000 





.0000 


0, 


,0000 


0.0000 


0.0000 


0.0000 


0.0000 


1.0000 


-0.0000 


-0.0000 


-0.0000 


-0, 


.0000 


CC15 


-0. 


.0000 


-0. 


.0000 


-0 


.0000 


-0 


.0000 


-0. 


0000 


0. 


.0000 





.0000 





.0000 





.0000 


0. 


,0000 


-0.0000 


0.0000 


0.0000 


-0.0000 


-0.0000 


1.0000 


0.0000 


0.0000 


-0, 


.0163 


CC16 


-0 


0000 





0000 





.0000 





.0000 


0. 


0000 


0. 


0000 





.0000 





.0000 





.0000 





0000 


-0.0000 


0.0000 


-0.0000 


-0.0000 


-0.0000 


0.0000 


1.0000 


-0.0000 


-0 


0379 


CC17 


-0 


.0000 


0. 


0000 





.0000 





.0000 


-0. 


0000 


0. 


0000 





.0000 





.0000 





.0000 


-0 


0000 


0.0000 


0.0000 


0.0000 


-0.0000 


-0.0000 


0.0000 


-0.0000 


1.0000 


0, 


.0000 


ES 


-0 


.0610 


-0 


0695 


-0 


.0507 


-0 


.0525 


-0. 


.0556 


-0. 


.0637 


-0 


.0648 


-0 


.0595 


-0 


.0430 


-0, 


,0495 


-0.0466 


-0.0487 


-0.0570 


0.0000 


-0.0000 


-0.0163 


-0.0379 


0.0000 


1 


.0000 



TABLE XXXIII: Statistical correlation coefficients from SNO's signal extraction for the salt phase night data. The numbering of the CC spectral bins follows the ordering of energy bins 
in Table XXX. 







NC 




CC1 




CC2 




CC3 




CC4 




CC5 




CC6 




CC7 




CC8 




CC9 


CC10 


ecu 


CC12 


CC13 


CC14 


CC15 


CC16 


CC17 




ES 


NC 


1. 


.0000 


-0. 


.3485 


-0 


.3437 


-0 


.3060 


-0. 


.2669 


-0. 


.2089 


-0 


.1571 


-0 


.1034 


-0 


.0545 


-0. 


.0317 


-0.0138 


-0.0067 


-0.0010 


-0.0008 


0.0001 


0.0000 


0.0000 


0.0000 


-0 


0693 


CC1 


-0 


.3485 


1 


.0000 





.1580 





.1400 


0. 


.1201 


0. 


0967 





.0671 





.0477 





.0237 





.0143 


0.0085 


0.0027 


0.0017 


0.0003 


0.0000 


-0.0000 


-0.0000 


-0.0000 


-0 


.0684 


CC2 


-0 


.3437 


0. 


.1580 


1 


.0000 





.1379 


0. 


.1183 


0. 


.0952 





.0661 





.0469 





.0234 


0. 


.0141 


0.0084 


0.0027 


0.0016 


0.0003 


0.0000 


-0.0000 


-0.0000 


-0.0000 


-0 


.0463 


CC3 


-0 


.3060 





.1400 





.1379 


1 


.0000 


0. 


.1048 


0. 


.0843 





.0586 





.0416 





.0207 


0. 


.0125 


0.0074 


0.0024 


0.0014 


0.0002 


0.0000 


-0.0000 


-0.0000 


-0.0000 


-0 


.0743 


CC4 


-0. 


.2669 





.1201 





.1183 





.1048 


1. 


.0000 


0. 


.0723 





.0505 





.0357 





.0178 


0. 


.0107 


0.0062 


0.0020 


0.0012 


0.0002 


0.0000 


-0.0000 


-0.0000 


-0.0000 


-0. 


.0527 


CC5 


-0. 


.2089 


0. 


.0967 





.0952 





.0843 


0. 


.0723 


1. 


0000 





.0404 





.0287 





.0143 


0. 


.0086 


0.0051 


0.0016 


0.0010 


0.0002 


0.0000 


-0.0000 


-0.0000 


-0.0000 


-0. 


.0483 


CC6 


-0. 


.1571 


0. 


0671 





.0661 





.0586 


0. 


.0505 


0. 


.0404 


1 


.0000 





.0199 





.0101 





0060 


0.0033 


0.0012 


0.0006 


0.0001 


-0.0000 


0.0000 


-0.0000 


0.0000 


-0. 


.0489 


CC7 


-0. 


.1034 





.0477 





.0469 





.0416 


0. 


.0357 


0. 


.0287 





.0199 


1 


.0000 





.0070 





.0042 


0.0025 


0.0008 


0.0005 


0.0001 


0.0000 


0.0000 


0.0000 


0.0000 


-0 


.0417 


CC8 


-0. 


.0545 


0. 


.0237 





.0234 





.0207 


0. 


.0178 


0. 


0143 





.0101 





.0070 


1 


.0000 





0021 


0.0012 


0.0004 


0.0002 


0.0000 


0.0000 


0.0000 


-0.0000 


0.0000 


-0. 


.0692 


CC9 


-0 


.0317 





0143 





.0141 





.0125 


0. 


0107 


0. 


0086 





.0060 





.0042 





.0021 


1 


.0000 


0.0007 


0.0002 


0.0001 


0.0000 


0.0000 


0.0000 


-0.0000 


0.0000 


-0. 


.0410 


CC10 


-0 


0138 





.0085 





.0084 





.0074 


0. 


.0062 


0. 


0051 





.0033 





.0025 





.0012 


0. 


.0007 


1.0000 


0.0001 


0.0001 


0.0000 


0.0000 


-0.0000 


-0.0000 


-0.0000 


0. 


.0016 


ecu 


-0. 


.0067 


0. 


.0027 





.0027 





.0024 





.0020 


0. 


0016 





.0012 





.0008 





.0004 





.0002 


0.0001 


1.0000 


0.0000 


0.0000 


-0.0000 


-0.0000 


-0.0000 


-0.0000 


-0 


.0359 


CC12 


-0. 


.0010 


0. 


.0017 





.0016 





.0014 


0. 


.0012 


0. 


0010 





.0006 





.0005 





.0002 





.0001 


0.0001 


0.0000 


1.0000 


0.0000 


0.0000 


-0.0000 


-0.0000 


-0.0000 





.0003 


CC13 


-0. 


.0008 


0. 


.0003 





.0003 





.0002 


0. 


.0002 


0. 


.0002 





.0001 





.0001 





,0000 


0. 


.0000 


0.0000 


0.0000 


0.0000 


1.0000 


0.0000 


0.0000 


-0.0000 


0.0000 


-0. 


.0229 


CC14 


0. 


.0001 


0. 


.0000 





.0000 





.0000 


0. 


.0000 


0. 


0000 


-0 


.0000 





.0000 





,0000 


0. 


.0000 


0.0000 


-0.0000 


0.0000 


0.0000 


1.0000 


-0.0000 


-0.0000 


-0.0000 


-0. 


.0139 


CC15 


0. 


.0000 


-0. 


0000 


-0 


.0000 


-0 


.0000 


-0. 


0000 


-0. 


0000 





,0000 





.0000 





.0000 


0. 


.0000 


-0.0000 


-0.0000 


-0.0000 


0.0000 


-0.0000 


1.0000 


-0.0000 


-0.0000 


-0. 


.0188 


CC16 





0000 


-0 


0000 


-0 


.0000 


-0 


.0000 


-0. 


0000 


-0. 


.0000 


-0 


.0000 





.0000 


-0 


.0000 


-0 


0000 


-0.0000 


-0.0000 


-0.0000 


-0.0000 


-0.0000 


-0.0000 


1.0000 


0.0000 


-0 


.0085 


CC17 


0. 


.0000 


-0 


0000 


-0 


.0000 


-0 


.0000 


-0. 


0000 


-0. 


0000 





.0000 





.0000 





.0000 





0000 


-0.0000 


-0.0000 


-0.0000 


0.0000 


-0.0000 


-0.0000 


0.0000 


1.0000 


-0. 


.0183 


ES 


-0 


.0693 


-0. 


.0684 


-0 


.0463 


-0 


.0743 


-0. 


.0527 


-0. 


.0483 


-0 


.0489 


-0 


.0417 


-0 


.0692 


-0. 


.0410 


0.0016 


-0.0359 


0.0003 


-0.0229 


-0.0139 


-0.0188 


-0.0085 


-0.0183 


1 


.0000 



TABLE XXXIV: Complete CC spectrum systematic uncertainties (in percent) from SNO's unconstrained fit signal extraction of the 391-day salt phase data set. The numbering of the CC 
spectral bins follows the ordering of energy bins in Table XXX. 



Uncertainty 


NC 


ES 


CCl 


CC2 


CC3 


CC4 


CC5 


CC6 


CC7 


CC8 


Energy scale (const.) 


3.8, 


-3.3 


1.9, 


-1.6 


-8.0, 


9.5 


-5.9, 


6.2 


-2.9, 


3.0 


0.6, 


0.0 


1.7, 


-1.8 


2.9, 


-3.9 


3.9, 


-3.6 


4.7, 


-5.0 


Energy scale (E dep.) 


0.1, 


-0.1 


0.1, 


-0.1 


-1.5, 


1.6 


-1.4, 


1.6 


-1.6, 


1.2 


-0.4, 


0.6 


-0.5, 


0.6 


0.0, 


-0.4 


0.9, 


-0.6 


1.3, 


-0.4 


Energy radial bias 


2.1, 


-2.0 


1.2, 


-1.1 


-5.7, 


6.1 


-4.3, 


4.2 


-2.1, 


2.0 


0.1, 


-0.1 


1.0, 


-1.2 


1.5, 


-2.5 


3.6, 


-3.9 


3.8, 


-2.9 


Energy resolution 


0.8, 


-0.8 


0.7, 


-0.7 


4.9, 


-4.9 


2.1, 


-2.1 


1.8, 


-1.8 


0.1, 


-0.1 


-1.2, 


1.2 


-1.8, 


1.8 


-3.9, 


3.9 


-2.6, 


2.6 


j8i4 mean 


-3.6, 


4.5 


1.3, 


-1.2 


7.1, 


-8.4 


7.0, 


-8.3 


6.3, 


-7.9 


4.4, 


-4.6 


3.3, 


-4.0 


2.3, 


-2.9 


2.6, 


-2.9 


1.1, 


-1.7 


/?i4 width 


0.0, 


0.0 


0.2, 


-0.2 


-0.9, 


-0.1 


-0.7, 


0.0 


-0.8, 


-0.3 


0.0, 


-0.2 


-0.2, 


-0.2 


0.0, 


-0.2 


-0.1, 


0.3 


0.3, 


-0.2 


Radial scale (const.) 


-3.0, 


3.3 


-2.6, 


3.0 


-2.6, 


2.5 


-3.4, 


2.3 


-1.7, 


2.4 


-2.4, 


2.6 


-2.5, 


2.3 


-2.8, 


2.4 


-2.6, 


2.9 


-2.5, 


2.8 


Radial scale (E dep.) 


-0.6, 


0.5 


-0.7, 


0.8 


0.2, 


-0.2 


-0.4, 


0.4 


-0.9, 


0.9 


-1.2, 


1.2 


-1.5, 


1.5 


-1.7, 


1.7 


-2.0, 


2.0 


-2.1, 


2.1 


Vertex x 


0.0, 


-0.0 


0.1, 


-0.1 


-0.3, 


0.0 


-0.2, 


-0.5 


-0.3, 


0.0 


-0.5, 


0.0 


-0.1, 


0.0 


-0.2, 


-0.1 


-0.1, 


-0.1 


0.1, 


-0.1 


Vertex y 


0.0, 


-0.1 


0.1, 


-0.1 


-0.2, 


-0.3 


0.0, 


-0.2 


0.1, 


-0.1 


-0.2, 


0.0 


0.1, 


0.0 


-0.3, 


0.0 


0.0, 


-0.3 


0.1, 


0.4 


Vertex z 


-0.2, 


0.2 


0.0, 


0.0 


-0.4, 


-0.2 


-0.9, 


-0.6 


0.4, 


0.1 


0.2, 


-0.6 


0.0, 


-0.3 


0.0, 


-0.1 


0.5, 


-0.2 


-0.3, 


-0.1 


Vertex resolution 


0.1, 


-0.1 


0.1, 


-0.1 


-0.4, 


0.4 


-0.2, 


0.2 


1.4, 


-1.4 


-0.3, 


0.3 


0.1, 


-0.1 


-0.6, 


0.6 


0.0, 


0.0 


0.4, 


-0.4 


Angular resolution 


-0.2, 


0.2 


5.1, 


-5.1 


-0.2, 


0.2 


0.4, 


-0.4 


0.2, 


-0.2 


0.1, 


-0.1 


-0.2, 


0.2 


0.0, 


0.0 


0.0, 


0.0 


-0.6, 


0.6 


Internal y 


0.0, 


-0.1 


0.0, 


-0.0 


0.4, 


-0.6 


0.3, 


-0.4 


0.2, 


-0.2 


0.1, 


-0.1 


0.0, 


0.0 


0.0, 


0.0 


0.0, 


0.0 


0.0, 


0.0 


Selection efficiency 


-0.1, 


0.2 


-0.2, 


0.2 


-0.2, 


0.2 


-0.1, 


0.2 


-0.1, 


0.2 


-0.1, 


0.2 


-0.1, 


0.2 


-0.1, 


0.2 


-0.1, 


0.2 


-0.1, 


0.2 


Backgrounds 


-1.0, 


0.0 


-0.0, 


0.0 


-8.5, 


0.0 


-0.1, 


0.0 


-0.1, 


0.0 


-0.1, 


0.0 


-0.1, 


0.0 


-0.1, 


0.0 


-0.1, 


0.0 


-0.1, 


0.0 




CC9 


CC10 


ecu 


CC12 


CC13 


CC14 


CC15 


CC16 


CC17 






Energy scale (const.) 


5.0, 


-4.2 


6.2, 


-6.7 


6.7, 


-5.8 


9.6, 


-7.4 


7.0, 


-6.7 


17.6, 


-9.7 


16.6, 


-17.4 


18.2, 


-10.2 


20.5, 


-17.5 






Energy scale (E dep.) 


0.7, 


-0.8 


1.3, 


-2.0 


2.3, 


-2.3 


3.9, 


-2.5 


2.5, 


-2.0 


9.3, 


-5.1 


9.1, • 


■11.9 


13.8, 


-8.8 


14.1, 


-12.2 






Energy radial bias 


2.8, 


-3.1 


4.1, 


-4.3 


4.7, 


-4.2 


7.0, 


-4.5 


5.7, 


-4.7 


12.4, 


-6.6 


11.8, 


-14.6 


13.5, 


-10.3 


13.2, 


-12.3 






Energy resolution 


-0.7, 


0.7 


0.0, 


0.0 


-2.1, 


2.1 


0.6, 


-0.6 


-2.3, 


2.3 


-5.8, 


5.8 


-14.8, 


14.8 


-9.9, 


9.9 


-17.2, 


17.2 






/3i4 mean 


0.6, 


-0.4 


0.2, 


0.4 


-0.1, 


-0.1 


1.4, 


-0.1 


0.3, 


0.1 


-1.0, 


4.4 


2.2, 


-1.4 


-5.9, 


9.0 


3.8, 


-2.1 






yS 14 width 


0.0, 


-0. 


-0.3, 


0.6 


0.0, 


0.2 


0.0, 


0.2 


-0.4, 


0.0 


-0.5, 


1.9 


0.5, 


-0.1 


-0.1, 


6.5 


-0.2, 


-0.2 






Radial scale (const.) 


-2.8, 


3.5 


-2.0, 


2.0 


-2.2, 


1.9 


-1.7, 


2.4 


-2.4, 


3.7 


-2.5, 


5.6 


-12.6, 


3.9 


-0.4, 


-1.4 


-2.8, 


2.9 






Radial scale (E dep.) 


-2.2, 


2.2 


-2.3, 


2.5 


-2.5, 


2.7 


-2.7, 


2.9 


-2.9, 


3.1 


-3.1, 


3.3 


-3.3, 


3.5 


-3.6, 


3.8 


-3.9, 


4.3 






Vertex x 


-0.1, 


-0.1 


0.0, 


-0.5 


-0.1, 


0.3 


-0.2, 


0.0 


0.7, 


-0.4 


0.2, 


-0.2 


-1.2, 


0.0 


-0.9, 


-0.6 


-0.2, 


0.5 






Vertex y 


0.0, 


0.1 


0.5, 


-0.2 


0.0, 


0.0 


0.5, 


0.0 


0.0, 


0.3 


0.5, 


-0.1 


-0.5, 


0.2 


-1.3, 


-0.4 


-0.1, 


0.1 






Vertex z 


0.7, 


0.0 


-0.2, 


-0.3 


0.1, 


-0.3 


0.7, 


0.1 


0.8, 


-1.0 


0.5, 


1.3 


-1.2, 


1.0 


-2.4, 


1.0 


0.1, 


-0.9 






Vertex resolution 


0.0, 


0.0 


-0.2, 


0.2 


0.1, 


-0.1 


-0.1, 


0.1 


1.0, 


-1.0 


3.3, 


-3.3 


-8.7, 


8.7 


-4.1, 


4.1 


0.8, 


-0.8 






Angular resolution 


0.0, 


0.0 


-0.6, 


0.6 


-0.7, 


0.7 


-0.8, 


0.8 


-0.8, 


0.8 


0.3, 


-0.3 


0.5, 


-0.5 


0.4, 


-0.4 


0.0, 


0.0 






Internal y 


0.0, 


0.0 


0.0, 


0.0 


0.0, 


0.0 


0.0, 


0.0 


0.0, 


0.0 


0.0, 


0.0 


0.0, 


0.0 


0.0, 


0.0 


0.0, 


0.0 






Selection efficiency 


-0.1, 


0.2 


-0.2, 


0.2 


-0.2, 


0.2 


-0.2, 


0.2 


-0.2, 


0.2 


-0.2, 


0.2 


-0.2, 


0.2 


-0.2, 


0.2 


-0.3, 


0.3 






Backgrounds 


-0.1, 


0.0 


-0.2, 


0.0 


-0.2, 


0.0 


-0.3, 


0.0 


-0.3, 


0.0 


-0.9, 


0.0 


-0.9, 


0.0 


-3.2, 


0.0 


-12.3, 


0.0 







